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Translocation of the Global Positioning System has proved to be a highly 
accurate method of position determination for onshore and airborne navi- 
gation, but it had not been previously evaluated for nearshore hydro- 
graphic surveys. The technique of translocation for hydrographic opera- 
tions involves the simultaneous reception of signals from the GPS 
Satellites by two independent receivers; one receiver onboard the survey 
vessel and one located at a known stationary position. A position 
correction is obtained at the stationary receiver, which is then applied 
to the shipboard receiver online or during postprocessing. An accuracy 
determination of translocated GPS was conducted at the Naval Postgraduate 
School, Monterey, CA, May, 1981. Two methods of positioning were used 
for comparison with GPS: 1) a least squares solution of three lines of 
position observed from three Wild T-2 theodolites, and 2) a position 
determined from ARTEMIS, a range-azimuth short range microwave position- 
ing system. Translocated GPS accuracies of 10 meters were determined. 
It is anticipated that greater accuracies will be obtained by using a 


more sophisticated receiver and more advanced processing methods. 
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Le  LHTEROLU CTION 


Shipboard operations for the purpose of oceanographic 
research, commercial transport, tactical operations, or 
hydrographic surveying require accurate and dependable 
positioning systems. The operation of most 2lectronic 
positioning systems is hampered by atmospheric effects on 
the propagated signal and design limitations which restrict 
operating ranges, positional accuracy and dependability. A 
universal positioning system that provides accurate 
information to the user regardless of atmospheric conditions 
and that is free of range limitations is desirable. ¥ith 
the advent of satellite positioning technigues, such a 
system is now available. 

The Navy Navigation Satellite System (NAVSAT), Or 
TRANSIT, has been operational since 1964 and has' proved 
dependable and sufficiently accurate for certain users. 
Although developed for military uses, the application of 
TRANSIT in the civilian community has continued to expand at 
an exponential rate. The system, however, is limited by the 
time interval between fixes, which is 90 minutes on the 
average, andthe requirement for precise knowledse of the 


ship's velocity. Although the system does not provide 





One method of improving system accuracy is by using two 
receivers Simultaneously in a technique terned 
translocation, or differential mode. The method is based on 
the assumption that signal propagation errors to two 
proximate receivers are nearly the same. One receiver is 
placed over a known position, and the variation of the 
observed position with time 1s determined relative to the 
absolute position. This variation is applied to the 
observed positions determined by a Simultaneously operating 
receiver located at another location--on a survey vessel, 
for example. The assumption underlying this technique is 
that the distance between the two units is not sufficient to 
introduce a significant additional propagation error. The 
resulting position of the translocated receiver should be 
more accurate than the position obtained using a ingle 
receiver. 

For hydrographic operations, the translocated Global 
Positioning Systen (GPS) provides an alternative to 
presently used positioning systems. The satellites orbit at 
20,000 kn, which theoretically would allow receiver 
separation of as much as 500 km without causing Significant 


additional system error. Weth this flezkibility, one 


ies 





sufficient positional accuracy for hydrographic surveying, 
it has been used for geodetic application by the Defense 
Mapping Agency (DMA) using fixed shore-based receivers to 
establish positioning networks in inaccessible areas 
worldwide. [{Ref. 1} 

The Global Positioning Svstem, or NAVSTAR, is a more 
sophisticated system designed to replace TRANSIT in the 
mid-1980s. Six satellites are currently in orbit with 
twelve more to be launched before the operational date of 
1987. The 18-satellite orbit configuration will provide a 
multi-user, passive system with 24-hour availability, global 
coverage, and possible fix intervals of less than one 
second. The operational network consists of three 
elements--the satellites, the ground-based control segnent, 
operated by the Space and Missile Systems Organization 
(SAMSO), and the user's receiver. Prototype recelvers have 
been tested under optimum conditions by SAMSO at the Yuma 
Proving Ground, Arizona, on both mobile and stationary 
platforms. Results show excellent system Stability and 
accuracy when compared to currently used positioning 


systems. (Ref. 2] 
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shorebased receiver could be placed at a central location 
and a survey of several hundred kilometers of coastline or 

ffshore area could be completed with only one additional 
receiver on the survey vessel. (Ref. 3] The value of this 
flexibility would be realized in the significant decrease in 
survey planning, cost, and time compared to that spent 
presently in conducting hydrographic operations. 

The accuracy of the system would be sufficient for all 
survey scales. The National Ocean Survey (NOS) defines the 
reguired accuracy of a positioning system to be no nore than 
one-third of the survey accuracy, which is 1.5 millimeters 
at the scale of the survey. {Ref. 4] Non-translocated GPS 
has been shown to be of sufficient accuracy for 1:80,000 
scale surveys and smaller. [Ref. 5] 

The object of our research is to determine whether GPS 
at the present operational level is adeguate for 1:10,000 
scale surveys. To do this, a position determined by two GPS 
receivers in the translocated mode is compared to the same 
position determined by two other methods. If translocated 
GPS can be shown to achieve specified inshore survey 
accuracies, the system would provide sufficient accuracies 


at all smaller scale surveys. 





Part of the difficulty in determining this level of 
accuracy concerns the reference system used. Most geodetic 
positions published by the National Geodetic Survey (NGS) 
are based on the North American Reference Datum of 1927 
(NAD-27) and the Clarke N66 wer slspsol dG , the best 
approximation to the shape of the earth for the continental 
United States. The satellite systems, however, use a 
reference system which is not based ona portion of the 
earth's surface. A mass-~centered ellipsoid referenced to 
the earth's center of mass has been computed using observed 
gravity and astronomical data. The best fit to this global 
ellipsoid is the World Geodetic System of 1972 (WGS-72), 
which is utilized by both “TRANSIT and GPS. Transformation 
equations exist which convert NAD-27 coordinates to WGS-72 
coordinates, but uncertainty in the gravity model produces 
errors as large as 10 meters due to the transformation 
precess alone. [{Ref. 6] 

To establish all positions on the same reference 
ellipsoid for this research, 2 TRANSIT receiver was acquired 
from DMA and a first order Doppler station was established 
at the central location of the test area. A third order 


geodetic survey was conducted from the TRANSIT station to 


a0 


establish two additional seal vons on the WGS-72 ellipsoid. 
Two Single channel Manpack GPS receivers were obtained from 
SAMSO. Stability testing and a receiver comparison test 
were conducted prior to shipboard testing tc determine a 
guantitative measure of receiver performance. One of the 
receivers was then placed over the TRANSIT station mark and 
the other aboard a survey vessel which maneuvered along 
predetermined track lines. Data were recorded via an 
interface to an #HP-9825 computer and printer at each site. 
Simultaneously, two additional independent measurements were 
made of ship position--one using a Trange-azinuth short 
range microwave positioning systen (ARTEMIS) and the other a 
position determined by a least squares solution of three 
lines of position observed from three 1-second theodolites 
set on the three pre-established WGS-72 locations. The 
comparison of the ARTEMIS and theodolite-derived ship 
positions with the translocated GPS ship position provided 


the test results. 
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II. EQUIPMENT DESCRIPTION AND OPERATIONAL CAPABILITIES 


Two positioning systems were chosen for comparison with 
the Global Positioning Systems a theodolite network anda 
microwave positioning systen. The theodolite network 
consisted of three shore-based units so as tn provide *three 
lines of position. The microwave positioning systen 
previded a range and an azimuth from a known location on 
shore. Both systems are accurate at short ranges and 
represent two contrasting methods of poSitioning--visual and 
electronic. The theodolite method of positioning at short 
range has been and continues to be used for hydrographic 
Operations as an alternative +o electronic systems when 
gecmetry and system accuracy limitations impose the need for 
a more flexible positioning systen. The requirement for 
experienced theodolite observers as well as the 
communication and logistics difficulties involved limit the 
use of this method to special cases. 

This chapter describes the ARTEMIS short-range 
electronic positioning systen, the operation of the 
sheodolite network as used during ‘the comparison test, and 


the TRANSIT satellite systen. 





> weet GoLsS POSITIONING SYSTEM 

The ARTEMIS position fixing system is a _ sShort-range 
\lcrowave systen (to 30 =kn), but “by Christiaan 
luygenslaboratorium B. V., in the Netherlands. It has been 
ised operationally (primarily in Europe) since 1972, when 
che system was introduced at the International Hydrographic 
conference in Monaco. At the present time, fifty systems 
ire in use. The system is unigue in design, is highly 
rccurate at short ranges, and is particularly useful for 
small harbor surveys, river surveys and relative positioning 
ises. Por this study, the ARTEMIS system, along with 
technical assistance, was obtained from the American 
iistributor, MARINAV Corporation of Houston, TEXAS. 

ARTEMIS is a range-aZimuth positioning system which 
2mploys a unigue tracking method. A stationary or fixed 
anit is placed over a known position on shore. Using a 
telescope that is mounted on the antenna, the antenna is 
Sighted onto another known point andthe known azimuth 
between the two locations is entered into the system via a 
digital display on the fixed unit. The unit is then locked 
into place for operation. A mobile unit is placed on the 


vessel whose position is desired. 
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The mobile and fixed antennas resemble radar antennas 
which, while operating, track each other so that the two 
antennas are always parallel to each other and perpendicular 
to the line of position between then. A microwave link is 
established between the two antennas during periods of 
measurement and information is relayed between the mobile 
and fixed units. (A voice link is also available during 
operation.) The fixed unit provides the direction and the 
azimuth of the mobile unit with respect to the point sighted 
on by the telescope. At the same time, the mobile unit 
transmits a coded interruption of the signal to the fixed 
unit, which replies with the same interruption. The tine 
difference between the transmission and reception at the 
mobile unit determines the distance between the two 
antennas. 

Distance accuracy is determined by the number of 
measurements averaged over time and is available in two 
modes--a dynamic mode which averages 1000 measurements or a 
static mode which measures 10,000 measurements ver displayed 
distance. The azimuth and distance information is displayed 
by the mobile unit so that personnel are not required at the 


stationary unit once operation has begun. The system can he 
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operated in long or short-range modes, which determine the 
amount of power output for the systen. 

The microwave frequencies used are 9.2 and 9.3 £4GHz. 
Problems arising with signal interference, observed with all 
microwave positioning systens, therefore can occur with 
ARTEMIS. The regions of range holes, Or miltipach 
interference, can be computed for up to 10 km range by prior 
knowledge of the antenna elevations above sea level and the 
operating distance from the shore-based antenna. [Ref. 7] 

Each antenna is designed with a 22° vertical beam width 
and a 2. horizontal beam width. A power source of 22 volts 
DC per unit is required, with an average current load of 2.5 
amps. 

The accuracy of the system depends on both the distance 
and the angle measured. The angle measurement has a 2 6 
probability of measuring within 2' of arc (0.033°). The 
distance measurement has a 2° probability that the mean 
distance value will vary within + 1.5 meters. 

One of the advantages of ARTEMIS over other electronic 
Systems is that the two lines of position defining the fix 
always intersect at 90. A unmique advantage is that the 


coverage around the fixed unit can be circular, with a 
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maximum range cf operation to 30 ka, dependent on antenna 
heights. The short range applications of the ARTEMIS fill a 
need not provided by Many currently used positioning 
systems. 

For this study, the ARTEMIS was used in conjunction with 
the positions derived fron the theodolit network +o 
determine the best estimate of the ship's position. The 
accuracy of the ARTEMIS was more than adequate for the 
ranges encountered during the test, and the ease of use and 
Operational versatility of the system made its use very 


advantageous. 


Be. THEODOLITE NETWORK 

A network of three Wild I-2 theodolites was used as the 
control positioning systen. The theodolite is widely used 
for surveying as a means of measuring angles between two 
observational locations--either ina horizcental plane for 
traverse or triangulation surveys, or a vertical plane for 
astronomical observations. 

The use of three theodolites aS a short-range 
positioning system was an accepted method of positioning 
ptior to development of alectronic distance measuring 


systems. There remain disadvantages when using electronic 
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positioning systems for hydrographic operations at short 
ranges (1-2 km). Most short-range electronic systems are 
designed for measurement of pesseron at =Pine-of=—sight 
ranges. This region includes a large portion of inshore 
hydrographic surveys. Inshore positioning involves nore 
thorough planning to obtain the required accuracy for the 
scale of the survey. For example, in a harbor or cove it 
may be difficult to place two range-range electronic units 
so that the area to be surveyed is not within the region of 
positioning system uncertainty, or 30° (radially) on each 
Side of the baseline between the two units. To reposition 
the electronic units as often as is required to obtain 
Optimum geometry is time-consuming and requires manpower and 
extensive planning, and can result in operational delays. 
Interference of the signal due to reflection from a calm sea 
surface, or multipathing, is also a problem at short ranges 
for microwave systems. 

Theodolites can be set at any intervisible location with 
preplanned geometry and can be used for the entire region to 
be surveyed. This methed requires less expensive equipment 
(three theodolites), more manpower (a recorder and an 


observer at each site), communication between the three 
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locations (via portable radio, for example) and.reasonable 
weather conditions. The basis of this procedure is that 
three angles observed From known locations to a 
predetermined point on the vessel constitute three lines of 
position. After removing all possible systematic errors and 
through a least squares solution of the three lines of 
position, the position of the vessel is determined. A major 
disadvantage is that observational errors are not determined 
in real time because the data cannot be processed 
immediately. 

The basic procedure for the method requires 
predetermined sites which are known relative to each other. 
This may be achieved by using published geographic positions 
or by conducting a local plane survey to establish new 
positions uSing preestablished stations nearby. Zither way, 
all three locations must be based on the same datum and 
preferably permanently positioned by a disk ora station 
mark. 

Each theodolite was aligned using a plumb bob or optical 
plummet over one of the known station marks or disks. By a 
predetermined scheme, each site was marked so as to be 


Clearly visible to the other sites (lights were used for 
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this study). The telescope of each theodolite was sighted 


onto one of the other sites and the horizontal circle set to 


read an arbitrary value close +o 0°. The crosshairs in the 
telescope were set to intersect the mark. When the 


theodolite was exactly sighted, the horizontal circle was 
read and recorded. (At this time the lights were turned 
Omi. ) The operation began when all three sites completed 
this step, called the "initialling procedure", or "initial 
pegntang' . The instrument was sighted on a point on the 
survey vessel that was visible to all operators. The ship 
was tracked by the operators by using the micrometer wheels 
on the instrument until a "mark" was called over the radio 
by the recorder at the control theodolite location. At that 
instant, the operator stopped tracking and the horizontal 
Circle was read and recorded. This observation constituted 
ae xX. After a series of fixes were recorded, or when the 
Survey vessel reached a predetermined location, a set of 
fixes was completed. The lights were then reactivated at 
each site, each observer intersected the original mark, the 
hMorsazontal citcle was again red@ andthe ‘final pointing" 
was recorded. If necessary, the instrument was relevelled. 
The process was then repeated by starting With the 


itnitialling procedure. 
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It was assumed that the theodolites were in good 
Operating condition. Collimation error had previously been 
determined. The observers were kept at the same locations 


throughout the test to minimize the error due to individual 


observer variations. Since only one portion of the 
herizontal circle was used, and no reverse readings 
observed, a systematic error was introduced. This was 


assumed to be minimized during computation of the vessel 
position. 

The observations were corrected for collimation error 
and for pointing error by evenly distributing the difference 
between the initial and final pointing values. The 
resulting angles were processed using an intersection 
computation with a least squares fit applied to optinize the 
result. The computations for this study were processed 
using an HP-9815 computer and software prepared by Cdr 
Ludvig Pfeifer, NOAA. The resulting position was used as 


the best estimate of the position of the ship. 


C. TRANSIT POSITIONING SYSTEN 
The TRANSIT, or Navy Navigation Satellite System, is the 
most widely utilized satellite positioning system at the 


present time. The network of satellites and ground tracking 
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stations has been in operation since 1964 and is becoming 
widely used by both the military and civilian communities. 
It is expected to continue operation until at least 1995, 
when the Global Positioning System is expected to replace 
ates 

The system consists of 2 series of TRANSIT satellites 
(ierg. 2.1.) in fixed polar orbits of 1095 km elevation, 
ferming a network within which the earth rotates (fig.2.2.). 
The orbits are spaced so that satellite passes occur, for a 
given ground location, at intervals of between 35 and 100 
minutes. Each satellite travels an observable distance of 
between 4400 and 7000 km per pass, providing a sufficient 
tine interval and baseline over which the user can record 
data. 

The idea of positioning an observer on the earth's 
surface by using the satellite as a reference systen 
Originated from observation of Sputnik 1, launched by the 
Soviet Union in 1957. By recording the change in frequency 
of the satellite signal as the satellite passed overhead, a 
Doppler shift could be measured. The Doppler shift gives a 
Spatial relationship as a function of time between the 


Observer and the satellite. Pf thee pesition of Eke 
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satellite is known at the time of observation, the 
observer's position can be determined. 

In the TRANSIT systen, the satellite positions are 
precisely determined by tracking stations located at 
Prospect Harbor, Maine; Rosemount, Minnesota; and Wahiawa, 
Hawaii. The tracking information (in the form of the 
Doppler frequency shift as a function of time) is relayed to 
the computing center at the Naval <Astronautics Group 
Headquarters, Po. . Srqu, Gamat Orniva. Here the actual 
satellite orbits are computed. This information is used to 
predict orbit information. Parameters of this oredicted 
orbit are relayed to two injection stations at Rosemount and 
Pt. Mugu. During the next satellite pass, the inforrmation is 
relayed to the satellite as part of the navigation nessage, 
With new updates occurring every 12 hours. 

The satellites therefore contain predicted orbit 
information which is received by the user via a satellit+t 
receiver. Several receivers are available with complexity 
and positional accuracy dependent upon the needs of the 
user. A computer is normally required to determine the 
user's mobile or fixed position in real time from a single 


Satellite pass. 
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For many users having a dual channel receiver and known 
velocity, this method of positicning provides sufficient 


accuracy with typical error of 27 to 37 meters root sun 


square. Optimum accuracy is obtained, however, if the 
receiver position is fixed, and several consecutive passes 
of the same satellite are recorded. During each pass, 


lasting from 10 to 16 minutes as the satellite travels from 
one horizon to the other, many Doppler counts (observations 
of shift in frequency as a function of tine) are recorded. 
After a predetermined number of passes are stored, the data 
can be forwarded to the Defense Mapping Agency in Brookmont, 
Maryland. There the precise ephemeris data from the ground 
tracking system is applied to the data from the recorded 
passes to provide the optimum solution. The resulting 
position has a typical error of 6.3 meters root méan square 
for a series of passes, with a 3-dimensional result of 1.5 
meters per axis repeatability after 25 precise ephemeris 
passes. This method is used to determine positions in 
remote areas by DMA geodetic survey groups using the 
Magnavox AN/PRR-14 Geoceiver (fig.2.3.). (Ref. 8] 

The TRANSIT system uses the World Geodetic System of 


1972 (WGS-7 2) reference ellipsoid which is based on the 
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Universal Space Rectangular Coordinate System--an earth, 
mass-centered coordinate systen. The wWGS-72 ellipsoid 
describes the best fit of the entire earth to a mathematical 
model. The satellite orbits are computed using this 
gecpotential model, the accuracy of which is based on 
gravitational and astronomical data. 

For the present study, a Magnavox AN/PRR-14 Geoceiver 
was acquired from the Department of Satellite Geophysics, 
Satellite Tracking Branch, DMA. The receiver was placed 
over a fixed point and 40 passes of one satellite were 
observed. The procedures described previously for an 
optimum solution of point positioning were followed. The 
resulting latitude, longitude, and height relative to the 
WGS-72 ellipsoid are coordinates of a first order TRANSIT 
position as defined by DMA standards. This position was 
used as the basis for horizontal control for the study (see 


test procedures). 


Sy 
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iii. Ges srolele DRSCRIPTION AND OPERATIONAL CAPABILITIES 











Reo ret) DESCRIPTION 


NAVSTAR, or the Global Positioning System, was designed 


for naxinun operational capability--high accuracy, 
twenty-four hceur access, passive user, and world-wide 
coverage. A brief description of the GPS, and the 


satellites (NAVSTARs), will serve as a background for the 
more in-depth discussion immediately following. 

GPS navigation is based on the measurement of four 
ranges, or pseudo-ranges, from each of four satellites with 
known, or predicted, positions. Each range determination is 
computed using the propagation velocity of the signal 
multiplied by the satellite-to-user travel time. Satellite 
transmission time is included as part of the navigation 
signal message. From the four known ranges, the user 
position and user clock error are computed on a three-axis, 
Orthogonal coordinate system. This computation is then used 
to transform the user position onto the WGS-72 ellipsoid. 

The remainder Of this Chapter will provide the 
background for GPS development and the three major systen 


segments of the GPS: the space system segment, *he control 
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segment, and the user system segment. Tncluded in the user 
system segment will be a discussion of the Manpack receiver 
as developed by Texas Instruments Corporation used for this 
study. 

With the successful operation of the TRANSIT system in 
the early 1960's, the Department of Defense (DOD) began to 
consider the specifications of, and the technology required 
for, the development and implementation of a second 
generation systen. One 17 April, 1973, the» Deputy Secretary 
of Defense stated that a Defense Navigation Satellite Systen 
was to be developed. {Ref. 9} Using the new synthetic 
oscillator technology, and the Department of Defense surface 
and nearsurface navigation objectives and requirements, a 
DOD interdepartmental task force developed the GPS as the 
second generation satellite navigation systen. With the 
United States Air Force (USAF) as the managing agency, the 
GPS became a joint service progran, With representatives 
from the DMA, Army, Navy, Marines, Coast Guard, and the 
North Atlantic Treaty Organization. [ Ref. 10] The Joint 
Programs Office was established within the Air Force Space 
and Missile Systems Organization to design, develop, and 


isplement the GPS progran. 





GPS was conceived and designed to provide continuous, 
real-time, world-wide navigation coverage for subsurface, 


surface and nearsurface operational vessels, vehicles, and 


amseratt. Accuracy of better than 10 meters was desired 
under all conditions, with navigational capabilities +o 


update and interface with other forms of existing and 
projected navigational systems. Another requirement of the 
GPS gave PGES Was nearly trouble-free and simple user 
operation, even in low Signal-to-noise environments. An 
encoded signal structure was required for the satellite 
system for security reasons. In addition, the signal 
transmitted by the NAVSTAR satellites had to include other 
operational information in conjunction with the navigational 
message, To provide for navigational requirements in event 
of hostilities, the GPS was designed to permit the use of 
redundant satellites. 

To achieve these objectives and requirements, the GPS 
program was organized according to three developmental 
tasks: concept validation, full-scale engineering 
development and systems tests, and production and operation. 
Starting in 1974, the concept validation task was initiated 


with the merger of the TIMATION series satellite program and 
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the USAF System 21B (a highly accurate three-dinensional 
Navigation system) progran. Tests with the Navigation 
Technology Satellite, NTS-1, utilizing two rubidium clecks 


(accuracy 5-10 parts per Tome ) and an orbit 7500 nautical 


miles (nn) above the surface and three ground based 
receivers (called inverted range) provided concept 
verification. {Ref. 11] Additional tests were conducted 


with an improved NTS-2 satellite positioned at an altitude 
of 10,980 nm and containing two cesium clocks {accuracies on 
the order of 1 part per nor) (Ref. 12] Orbital height of 
10,980 nm permitted the evaluation of the effect of the 
planned GPS constellation altitude on navigational 
accuracies using the improved cesium frequency standards. 
In 1980, NTS-3 was launched and the Advanced Developmental 
Model (ADM) hydrogen maser frequency standard evaluated for 
use on the GPS (providing an expected frequency accuracy of 
1 part per 10’°). [Ref. 13] The results obtained from the 
NTS-3 satellite have not yet been released. Since 1978, 
additional satellites have been orbited to provide the 
Six-satellite constellation required to fully evaluate the 


first development task--concept validation, Phase I. 
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During Phase I, user requirements were consolidated to 
allow the development and construction of four basic types 
of test receivers: continuous tracicsing receivers, 
sequential tracking receivers, lcw cost receivers, and the 
Manpack receivers. (Ref. 149] These receiver types were 
then constructed for Phase I testing by various contractors 
and their performances evaluated aboard helicopters, 
aacreraft, Ships, vehicles and personnel at the Joint 
Programs Office's testing facilities. Phase I, concept 
validation, was completed in late 1979 and Phase II testing 
was begun in early 1980 with full scale engineering, 
developmental, and systems tests. 

Phase II tasks, finishing the operational and control 
segments of the testing, are scheduled for completion in 
mid-1983. The operational control station will be relocated 
at aocentral continental United States site during this 
phase. Additionally, all satellite monitor receivers at the 
remote monitoring stations will be updated to operational 
status. User equipment and receiver selection will be made 
in two increments. Pour contractors will be asked to 
compete for the second increment of the receiver selection 


process by providing upgraded and refined user receivers for 
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extensive testing and evaluation. m7) of tie LOM 
weracractors Will then be selected to continue the 
competition and will provide additional receiver design, 
refinement and modificaticn, leading to prototype production 
and extensive field testing. 

When Phase II is completed in mid-1983, Phase III, 
production, will commence with the selection of one of the 
two contractors to produce the user receivers. Manufacture 
of these receivers will commence shortly after selection, 
and commercial availability of recaivers is expected in 
1985. {Ref. 15] The GPS should be fully operational with 
world-wide coverage using a constellation containing a total 
of eighteen satellites in 1988. [Ref. 16, 17] 

The GPS is considered to be composed of Ghee 
inter-related system segments: the space system segment, 
the control segment, and the user system segment. Of these, 
the space system segment may be considered the nost 
important since the other two require the space segment to 
be operational at all tines. Operationally, the space 
system segment consists of 18 satellites placed in orbits at 
altitudes of 10,898 nm and having orbital periods of about 


12 hours. [({Ref. 18] 
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The original plan for the satellite constellation 
consisted of 24 satellites in three planes of 8 satellites 
each. These three planes were to be placed 45° of longitude 
apart and at an inclination of 63° to the equator. 
(Ref. 19) As aresult of 1980 budget cutbacks the space 
system segment has been reduced to a total of 18 satellites. 
[ Ref. 20] Several different types of 18 satellite 
constellations have been proposed, but all have drawbacks. 
One which appears the most promising is the 'Nonunifornm 18° 
constellation containing 6 planes, each of which contains 3 
satellites spaced 120 apart (fig.3.1.). The areas of user 
outages (fig.3.2.) are a disadvantage of this configuration. 
The user outages are indicated by the darker areas shown on 
Pigure 3.2. These outage areas will appear to retate about 
the globe at the indicated latitudes, and positioning with 
the GPS will be intermittent, or of decreased accuracy for 
users at these latitudes. ane finiad constellation 
configuration is to be selected early in 1982. Ref. 21] 

Each NAVSTAR Satellite has the following 
characteristics: weight of 982 pounds, length of 17.5 feet 


From tip to tip of the 5 square solar panels, ‘three atomic 
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ia. 3.1. The Nonuniform 18 Sateliite Consteliation 
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frequency standards (clocks), and three nickel-cadmiun 
petocries for operations in darkness and during peak load 
periods. Stabilization of each NAVSTAR is maintained by 
four skewed reaction wheels which produce three-axis 
stabilization. NAVSTAR station-keeping and momentum damping 
capability is provided by a hydrazine propulsion systen 
Within each satellite. Life expectancy of each unit is five 
years, and expendables are expected to last for up to 7 
years with normal usage. 

Pigure 3.3. illustrates the parts of each NAVSTAR. A 
twelve element shaped-beam helix antenna transmits the 
navigation message on two L band frequencies to user systen 
segments. Simultaneously it receives satellite status, 
clock corrections, ionospheric data, and ephemeris constants 
from the control segment on two S band frequencies. 
(Ref. 22] The two coherent L band frequencies are centered 
on 1227.6 MHZ (L1), and 1575.48 MHz (L2). (Ref. 23] 
{Ref. 24] Both frequencies are necessary to determine the 
total ionospheric delay correction for Signal propagation 
time (Table I). Errors produced by the lonospheric delay 


in signal velocity contribute most of the error in the 


pseudo-range. BY utilizing the Lil and L2 frequencies and 
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Fig. 3.3. NAVSTAR Satellite 
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TABLE I. Ionospheric Group Delay as a Function of Frequency 


GD 
(neglecting third order and higher terms which contribute minimal 


error) 
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the following algorithm, the error is decreased but not 


totally eliminated. [Ref. 25] 


= Tent aa. 


Where: _ 
E GDL, Be 
Ly 


where A = ionospheric constant and * GOL, ¢ T GDL, are the 
respective ionospheric group delays at frequencies Li and 
L2. 80th frequencies are also modulated to transmit the GPS 
time ccdes and the navigation message. 

GPS time is transmitted using two pseudo-random noise 
chip codes--the Coarse Acquisition (C/A) and the Precision 
(P) codes. Modulation on the L1 frequency carries both the 
C/A and P codes in phase quadrature, while the L2 frequency 
Carries only the P code. These codes identify each acquired 
Satellite by matching the unique pseudo-random noise code 


pattern generated by each satellite with similar user 


generated codes. 
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Both codes are required to measure the navigation signal 
propagation time. The phase shift is measured and used to 
match the user generated signal with the incoming satellite 
Signal. (Ref. 26} The C/A code has a repetition period of 
one millisecond and has a short code stream, transmitting 
1.023 million bits per second (bps). This code is easily 
acquired and matched with the user generated pseudo-randon 
noise code but only provides a coarse time signal for 
computing the pseudo-range. However, the P code has a long 
code streaan, repeating only every seven days, and it is 
transmitted at the higher rate of 10.23 million bps. With 
this higher data transmission rate and iong repetition 
period, a more accurate time signal is available for 
pseudo-range determination. However, the long period and 
high data rate makes the P code difficult to match with the 
user-generated pseudo-random noise code pattern unless the 
receiver has a highly accurate time standard, and the 
approximate receiver location is known. In litem of this, 
the usual method of acquiring the P code is by matching the 
C/A code of the desired NAVSTAR, then using the Handover 
Word (HOW), to locate the correct P code seguence within the 


7-day pseudo-random noise code. GPS time is transmitted 
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with the HOW every six seconds and provides very accurate 
GPS time to the user. Both of these codes provide the 
proper time reference for computation of the pseudo-range. 
The C/A code provides a less accurate time reference than 
the P code. The selection of the code depends only on the 
navigational accuracy desired. 

The navigation message is also modulated on frequencies 
L1 and L2 but ata much slower rate (only 50 bps) than 
either the C/A or P_ code. On both frequencies, the 
navigation message is identical and consists of 5 subsets of 
6 seconds duration each. Each navigation message begins 
with the telemetry message (TLM) and the HOW message, which 
are used to transfer from the C/A to P codes. The remainder 
of the entire navigation message allocated to each 
subsection is transmitted after the TLM and HOW messages. 
Information in the navigation message includes approximate 
satellite ephemerides, status of other NAVSTAR satellites, 
parameters for clock corrections, atmospheric signal delay 
corrections, satellite performance status, momentum dump 
status since last upload, time since last upload, almanac 
and identification codes for all satellites, and provisions 


for inclusion of any special or important messages from the 





control system segment to the user. [ Ref. 27, 28] 
Navigation is thus achievable by combining the information 
transmitted by the space system segment in the form of the 
C/A and P codes and the navigation message. 

Accuracy and proper function of GPS are governed by the 
second segment--the control segment. The control segment 
monitors and uploads the individual NAVSTAR satellites, 
checks operational characteristics, determines the ephemeris 
and almanac data and corrections required, computes the tine 
delay and clock corrections for each satellite atomic 
standard, determines the atmospheric delay corrections for 
pseudo-range computations and includes special messages in 
the body of the navigational message. These functions are 
under the command and operational control cf the Master 
Sentrol Staticn, currently located at Vandenberg AFB, 
@altfornia. Four unmanned remote monitor stations, located 
in Hawaii, Guam, Elmendorf AFB, Alaska, and Vandenberg AFB, 
collect the pseudo-ranges to each satellite, the change in 
pseudo-range of each satellite signal, local meteorological 
data, and the remote monitor station atomic standard 
parameters. Each set of data is transmitted, on request, to 


the Master Control Station, where all the data is combined 
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and the satellite upload corrections and parameters are 
computed. These corrections and parameters are then relayed 
te the Upload Station, Vandenberg AFB, for transmittal to 
the respective NAVSTAR satellites. The uploads occur when 
the satellites are first visible to the Upload Station and 
when the User Range Error exceeds 4 meters, as computed by 


the Master Control Station. (Ref. 29, 30, 31] 


The last segment, the user system segment, combines the 
information provided to the space system segment (by the 
control segment) with the information received from the 
Space system segment to obtain the pseudo-ranges from four 
NAVSTAR satellites. These ranges are computed using the 


following algorithms: 


a = ena t,. = ake 
where R,; = true slant range to satellite i, c = speed of 
light in a vacuum, tp = time of received GPS signal by user 


receiver (assumed to be simultaneous from all satellites), 

= tine of transmission at satellite ia and 
i 

St, = atmospheric propagation delay time from satellite i 
] 
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(fig.3.4.). Assuming a user clock error, the true range can 


be expressed in terms of a pseudo-range: 


P = R. + —_" + — 
where a = pseudo-range including clock error, At, = time 
delay in the user receiver cf the received signal, and 
aS: = time delay in satellite transmission from satellite 


We 


To determine the slant ranges: 


Z 
By = Mes -2) ‘ (Ys -1)° #(z_ -2)° 


; 


R = xs =n ts n° + (25 -2)° tcht, +c( dt, -At. ) 


1 1 


The following quantities are known: yey, Zz obtained 


from navigation ephemeris data; A ty computed using the 


1 
icnospheric delay correction Obtained fron observing 


frequencies L1 and L2 tropospheric delay neglected); A a: 
obtained from the navigation message as the satellite clock 
SEror. Thus, there are four unknowns: xX, JY, z (user 
position coordinates) and 4Atp ( user clock error or bias). 
To solve this set of equations, four satellites must be 


observed simultaneously. [Ref. 32] 


55 





Should only three satellites be observable, the altitude 
may be assumed and ene solutions can be computed for the 
other two location variables and the user clock error. 
Sclutions of this type (altitude hold) significantly degrade 
position accuracy. 

The above description is over-simplified since the 
relativistic effects of NAVSTAR satellite velocity, as 
observed by the user, are neglected in terms of the 
satellite time delay correction and the atmospheric 
propagation delay correction. The actual user velocity, in 
terms of velocity inthe x, Y, sad Zz dine@tions, is 
calculated using the Doppler shifts of the two frequencies 
L1 and L2 and the estimated satellite position taken fron 
the ephemeris data. With the relativistic effects as an 
integral part of the data set, solution of the navigation 
position problem becomes more complex. Phase I receiver 
manufacturers have included them in the solution of the 
problem by several different methods. Methods of solution 
depend on the usage of the receiver (whether static, high or 
low dynamic) and the expected cost/benefit ratio for the 
receiver (in terms of the total cost, weight and size). 


Current state-of-the-art design technology and navigational 
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algorithms, developed since Phase I, will be included in the 
Phase II receivers and will significantly improve accuracy 


of future receivers. 


B. MANPACK RECEIVER 

The Phase I receiver used for this research was the 
Texas Instruments GPS Manpack, also known as the Manpack 
Vehicular User Equipment (MVUE). MVUE is a single channel 
low dynamic receiver, usable for receiver velocities from 0 
to 25 meters per second. [BRef. 33] System configuration is 
illustrated by figure 3.5. Composing the system are the 
Manpack control data unit (CD00), antenna and vehicle nount, 
and vehicle power converter. The Manpack is the receiver 
meee (229. 3.6.) with provisions to connect an external data 
output, through the MVUE Instrumentation System (MIS) 
interface, figure 3.7. Pacilities for an internal 
preamplified antenna nount, or an external preamplified 
antenna, and an attachment for either an external power 
source or the integral battery pack are also shown (figure 
3.7). Together, the MVUE and the S$BP/9900 microprocessor 
data processing unit have storage capacity for 47,360 
Sixteen bat words. This storage Capacity allows 


interrogation of NAVSTAR satellites sequentially, with the 
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receiver computing a position fix every 24 seconds. The 
relativistic velocity effects on the received signals are 
determined by matching the incoming frequencies with the 
expected frequencies. The user velocity relative to the 
Satellite is determined from the cbserved frequency shifts. 
The nominal frequency shift for each satellite transmission 
is computed using the simultaneously transmitted ephemeris 
data. Total weight of this receiver unit is 33.5 pounds, 
and it has a volume of 17,958 cubic centimeters. 

The second element of the MVUE, the CDU (fig. 3.8.), is 
the input/output device for the non-automated Manpack user. 
Inputs are entered by keying the pressure pads set in six 
rows of four pads per row. To obtain a navigational 
position, the approximate user latitude, longitude, and 
Greenwich Mean Time (GMT) are keyed. After the warmup 
period of 3 to 5 minutes, the satellite signals are 
received, and the CDU will display the actual position. 
Additional input may include NAVSTAR satellite selection, 
user altitude, and other position datums. Output options 
include position, GMT ‘time, and system error messages. 
Other inputs and outputs are listed in the MYJE manual. 


[ Ref. 34] 
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Automated users may elect output either to a radio 
interface or through the MIS interface. For this research, 
a special interface unit, designed and constructed at the 
Office of Testing and Evaluation of Systems, NOAA, was used 
to link the MIS interface board within the receiver with a 
Hewlett Packard 9825T computer. With this interface, user 
latitude, longitude and GPS time were logged onto magnetic 
tape for future processing. Figure 3.9. and Table II are 
the engineering documents used in construction of the MVUE 
interface unit. 

Other MYUE elements are the antenna andthe antenna 
external preamplifier vehicle mount. Located within the 
antenna are two antenna Subassemblies associated with 
frequencies L1and L2. Both of the subassemblies are 
connected to the external preamplifier with two lengths of 
coaxial cable, enclosed in a flexible conduit. The entire 
unit, antenna and preamplifier, is connected to the Manpack 
receiver with a power cable and a coaxial antenna cable. 
Cable length is determined by the location of the external 
antenna but is limited by the minimum strength reguired by 
the Manpack receiver. The combined antenna and preamplifier 


unit is 36 inches long, and it weighs 1.9 pounds. 
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Wire Color 
Code 


90 
9] 
o2 
o3 
94 
5 
96 
2). 
934 
e230 
936 
oGy 
945 
946 
947 
948 
98 
927 
928 
901 
SIZ 


TABLE I’. MANPACK Interface Wiring Code 
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The power converter must suoply 2% volts DC to the 
Manpack receiver. During this research, a power supply 
which converted 120 volts AC to 24 volts DC was used. This 
power supply was considered more stable than the power 
converter provided with the MVUE, and the voltage output and 
power were controlled and monitored during the testing 


period. 


C. DIFFERENTIAL GPS 

As mentioned earlier in this chapter, atmospheric 
propagation delays are the major source of error in the 
Navigational equations because ionospheric propagation is 
poorly understood. The tropospheric propagation delays are 
neglected in the navigation equations, and the resultant 
error is dependent on local temperature, pressure, humidity, 
and NAVSTAR elevation and azimuth. TIropospheric corrections 
are usually neglected if the observed satellite is more than 
5° above the local horizon. [Ref. 35] 

Additional error sources are receiver and satellite 
clock perturbations, lengths of signal paths from different 
Satellites, multipathing delays, satellite ephemeris errors, 
and receiver noise. Combined, these factors lead to errors 


of 3.6 to 6.3 meters, as table III illustrates. [Reft;: 36 ) 
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TABLE III.RANGE ERROR BUDGET 


Uncorrected Error Source User Equivalent Range Error, 
Cis) 


Feet meters 
SV Clock Errors = 
ephemeris Errors a0 lias 
Atmospheric Delays SeU iy 20 C= Sec 
Group Delay (SV Equipment) ere 126 
Multipath 4.0 - 9.0 eer) 
Receiver Noise and Resolution 50 1.5 
Vehicle Dynamics " 
Root Sum Square Hi eoe=ecOe/ BO = OS 


*Two nours after update 





One method used to reduce the magnitude of this error is 


that cf differential, or translocated, positioning. To use 
this method, one receiver is positioned over aie known 
geographic position. Pseudo-ranges to each satellite are 


determined, and a differential correcticn is obtained for 
each pseudo-range. These corrections are then applied to 
the satellite pseudo-range observed at a second mobile 
receiver. With the differential corrections applied, the 
net error of the mobile GPS receiver is significantly 


reduced. [{Ref. 37] 
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iV. 23ST PROC EDeRES 


A. POSITIONING NETWORK 

One of the considerations in planning this test was to 
examine the various methods of positioning the shore-based 
equipment. It was desirable to use closely spaced (1-2 km) 
locations within the immediate surroundings of Monterey 
harbor for reasons of ready access by personnel during 
installation of equipment and for optimum accuracy of the 
positioning systems chosen for control. AS mentioned 
previously, the positioning systems chosen for the 
comparison (the three theodolites and the ARTEMIS) give best 
results at relatively short ranges. 

Two predetermined station marks were selected for the 
test--1) USE MON 1978, a bronze disk set by the Army Corps 
of Engineers, located on a sand dune on the Naval 
Postgraduate School property 400 meters from the ocean at an 
elevation of 40 meters above sea level and 2) MUSSEL 1932, a 
bronze disk set by the Coast and Geodetic Survey in a rock 
outcrop 7 meters above sea level at the shoreline on the 
property of the Hopkins Marine Laboratory of Stanford 


University. The third station mark was set in the center of 
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the Monterey Coast Guard Pier for purposes of the test. 
This was a DMA bronze disk stamped GEOCEIVER STA 31370. 
This location was chosen for access to power facilities 
available at the pier and because of central location to the 
test area. 

The position of the Geoceiver station was established 
using a Magnavox Geoceiver (Model AN/PRR 14) obtained from 
the Department of Satellite Geophysics, Satellite Tracking 
Branch, DMA/Hydrographic Topographic Center. The equipment 
was operated in accordance with DMA requirements to obtain a 
documented first order station. The data were submitted to, 
and processed by, DMA/HTC in Brooknmont, Maryland, using 
precise ephemeris information to provide the WGS-72 
poSition. Michael Ellett of DMA provided the technical 
expertise for acquisition of the data and maintenance of the 
gecceiver during the test period. 

A third order, Class I, closed traverse was conducted 
following the standards of the National Geodetic Survey 
(NGS), using a Wild T-2 theodolite to measure angles and a 
Motorola MRA V Tellurometer to measure distances. The 
TRANSIT position computed by DMA/HTC was used for the 


initial position of the survey, anda Polaris observation 


/| 





using the theodolite provided the initial azimuth (fig. 
Beet.) . Reetc@tat CLOSUre GErEOr in distance of 0.061 feters 
and angle of 5.4" was achieved. The positions determined 


from the traverse and used for the control network are as 


follows: 

USE MON 36 36°08 .424"N 121 52°40.301"8 
MUSSEL 36 37°21.882"N 127 54'16.048"W 
GEOCEIVER STA 31370 36 36'36.246"N 127 53*29.693"8 


B. POSITIONING EQUIPMENT 

AS previously described, three positioning systems were 
used during the test--the three theodolites, the GPS, and 
the ARTEMIS. In addition, the geoceiver was set to acquire 
data during the shipboard portion of the test. During the 
time of day that the GPS satellites were available (Table 
IV) only one or two observations per day of the same 
TRANSIT satellite were possible. The geoceiver data were 


recorded but were not used for this study. 
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1. Theodolites 





This section describes the placement and operation 
of the equipment during the test. One theodolite was set on 
a — at each of the three stations, and each was 
occupied by an observer and a recorder. Communication 
between stations was maintained via portable radios. 

Because tests could be conducted only between 10 PM 
and 2 AM local mean time (LMT) due +0 GPS availability 
times, it was necessary to plan a systems of lights for the 
stations. The three theodolite sites were made visible by 
placing battery powered lights directly over the station 
marks. In addition, battery operated light packs for the 
theodolites were used. All shorebased lights and light 
packs for the test were provided by the Pacific Marine 
Center, NOAA. 

The ship system was designed with an omnidirectional 
pair of lights anchored around the support pipe for the GPS 
antenna. Also, the equipment shelter on the boat deck of 
the ship was lighted by spotlights so that the GPS antenna 
#as easily seen by the shorebased theodolite operators. 

During the test, each theodolite observer initialled 
on one of the other stations and recorded the reading on a 


leg sheet. The theodolite operators then intersected the 


is 


meghts One the ship GPS antenna using the theodolite 
telescope and tracked the ship until a 'mark' was called 
over the radio from the personnel aboard the ship. The 
angle measured at this time was recorded on the log sheet as 
Xs The time interval between fixes was approximately 
one minute. Each series of approximately twenty fixes 
constituted a set. At the end of each set, a final reading 
was made of the angle to the reference station. The fixes 
were numbered consecutively throughout the test, and the 
Greenwich Mean Time (GMT) of each fix was recorded on the 
ship log. sheets. The length of the sets depended on the 
ship position and heading. 

To correct the observed theodolite angles, the mean 
et the final pointing «nd the danitdiad pointing was 
determined for each set and for each station. This value 
Was assumed to be aoconstant throughout the set and was 
removed from each angle. A collimation error for each 
theodclite was also removed from each theodolite 
observation. The same observers and theodolites were used 
nightly at all stations except MUSSEL, where two different 
individuals observed over the test week. This continuity of 
personnel and equipment was designed to reduce operator and 


instrument error to a mininmun. 
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2. eebS and ARTENTsS 
To colocate the different pieces of equipment at the 


GEOCEIVER STA 31370, a wooden tower was built as shown in 


figure 4.2. This tower was set over the station for the 
period of the test. An NPS equipment vehicle was parked as 
close as possible to the tower. Electrical power was 


obtained from the Monterey Coast Guard facility. 

a Prior to the ship test, both GPS antennas were 
placed on the tower and centered over the station mark. For 
eye first two nights, 9 and 10 May, equipment problems 
allowed the use of only one system each night and data were 
acquired only to determine the stability of each systen. 
During the next night, 11 May, position outputs for the two 
Simultaneously operating GPS receivers were recorded. Tis 
Static test gave a determination of any difference between 
the two receivers under almost identical conditions. 

For the ship test, the GPS antenna was centered over 
the mark and was anchored at the top of the tower; the 
geoceiver antenna was centered over the mark and anchored at 
the middle level; and the theodolite on the tripod was set 
over the mark at ground level. This arrangement provided 


colocation of all three systems. It was assumed that no 
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error waS introduced by the presence of the wooden 
secucture. The shorebased ARTEMIS antenna unit was offset 
from the theodolite to eliminate interference in the signal 
due to the tower. An offset position was computed relative 
to the station mark and was used in postprocessing of the 
ARTEMIS data. 

Except for the antennas for the GPS, the ARTEMIS, 
and the geoceiver systems, all other equipment was installed 
in the equipment vehicle, where an operator was stationed 
throughout the test to monitor the data quality. 

The GPS Manpack is a manually operated systen. As 
mentioned previously, it was decided that manual operation 
of GPS was detrimental to the goal of the test, and that an 
interface that would provide an automated printout of 
positions from the Manpack as well as a printout of GPS time 
would be desirable. Knute Berstis and Gary Whitsell of the 
Office of Testing and Evaluation of Systems (OTES), NOAA, 
assembled two interfaces that provided data storage on 
cassette tapes as well as the desired printout. Each system 
Operated in conjunction w#ith a HP-9825 computer and a 
HP-7845 printer/plotter. Software was provided by OTES, 


NOAA. This arrangement permitted storage of the satellite 


19 





messages for further processing and provided an automated 
printout of the ARTEMIS data. A GPS Manpack had been sent 
mo OTES for testing during the construction of the 
interface. Technical expertise and Operation of the 
interface during the test were provided by OTES. 

The data from the GPS receiver were printed at an 
interval of 24 seconds, and consecutive data were 
continually recorded on tape and printed during the test 
periods. In addition, a manual record of the GPS position 
from the CDU was logged as the 'marks' were called via the 
radios. 

The ship positioning systems were installed as shown 
in figures 4.3. and 4.4. An equipment shelter was secured 
to the boat deck of the R/V ACANTA. A Mength of six-inch 
diameter polyvinyl chloride (PVC) pipe was attached to the 
aft wall of the building on the centerline of the ship. The 
GPS antenna andthe system of lights previously nentioned 
were set on top of the pipe. 

The ARTEMIS mobile antenna unit was set on a 
platform located on the centerline of the ship and at the 
top of the building high enough to not interfere with the 


Signal from the ship radar bean. For the vsveriod of the 
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test, the weather conditions allowed Ship operations without 
use of the radar, thus there was no possibility of damage to 
the ARTEMIS unit. 

All other equipment--the ARTEMIS receiver unit, the 
GPS receiver, control data unit, power filter, and 
interface, the HP-9825, the HP-7845 plottery/ printer, and 
the real time clock--was placed inside the building. An 
Operator monitored the equipment during the ship tests and 
another individual radioed the ‘'mark', Or Birt, to the 
theodolite operators on shore. The equipment used at the 
shore location was installed ina vehicle and consisted of 
the same components as on the ship. Operation of equipment 
during the ship tests was monitored by an individual located 


in the vehicle. 


C. LOG OF EVENTS 

The test was divided into two portions--1) a shorebased 
test which consisted of two nights of stability testing, and 
one night of receiver comparison tests, and 2) the ship test 
consisting of three nights of shipboard operation. The R/V 
ACANIA from the Naval Postgraduate School was the platfora 


used for the shipboard tests. 
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ie Stale hisy Test 
Operations began on 9 May, 1981, GMT. One of the 
GPS systems was placed over the station mark, and continuous 
data was recorded during the satellite availability periods. 
In this manner, system 1 was tested on 9 May, and system 2 
was tested on 10 May. A position was defined as a quality 
observation if four or more satellites were used by the 
receiver to obtain the output. 
mee 0oCcotatic Test 
Both GPS antennas were centered over GEOCEIVER STA 
51370, and data from Simultaneous operation of both 
receivers were recorded at the equipment vehicle. One 
antenna was placed on top of the wooden structure and the 
other immediately below it on the middle level. The 
vertical displacement was assumed to cause negligible error 
in the horizontal positions. Data were acquired on 11 May 
frem both receivers with a quality observation defined as a 
position determined from four or more satellites. 
3. Shipboard test | 
Por the entire period of shipboard testing, the wind 
speeds were less than 5 knots and the seas were calm with 1- 


2 feet or less of swell. Visibility was excellent throughout 
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the test week. The ship track for each set is shown in 


figure 4.5. A log of subsequent operations follows: 


12 May---The GPS Manpack receiver was installed on the R/V 
ACANIA and GPS and theodolite operations were conducted. 
The ship was allowed to drift during sets 1-3 (fix 1-69) 


With 23 fixes per set, in a east-northeasterly direction. 


13 May---The ship cruised at 4 knots on track lines from 
south to north for sets 4 and 6, and lines from north to 
south for sets 5 and 7. During fixes 139 through 157 the 
radio at the Coast Guard pier was disabled, and no data were 
acquired by either the theodolite observer or the operator 


in the equipment vehicle. 


14 May---The shipboard GPS system was inoperable. The 
trouble was traced to a faulty antenna cable, requiring 


replacement. No data were acquired. 


15 May---The ARTEMIS, the GPS andthe theodolite network 
were allused for this test period. Data from sets 8 


through 13 were acquired (fixes 165-262). 


16 May---The control data unit on the shore based GPS systen 
failed before the start of ship operation and could not be 


repaired without factory maintenance. 
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A total of 252 fixes were recorded by the theodolite 
network and all GPS positions corresponding to the 
theodolite fixes were manually logged. The ARTEMIS azinuth 
and distance information and consecutive GPS positions on 
the ship and in the vehicle were stored on cassette tapes 


and printed by the HP-7845 printer/plotter. 





Ve. DATA PROCESSING METHODS 





Processing of the thesis data was divided into four 
interrelated, but separate, tasks. These tasks were the 
basic geodetic survey, theodolite and ARTEMIS ship position 
determination, differential GPS ship position determination, 
and the statistical analysis of the ship position 


ccmparisons. Each task will be discussed separately. 


A. BASIC SURVEY COMPUTATIONS 

As described earlier in the Test Procedures, a basic 
geodetic Survey of the test area was conducted +o place the 
theodolite control stations on the World Geodetic System of 
1972, or WGS-72. To begin the survey, a Magnavox AN/PRR 14 
Geoceiver was placed over the center survey disk, GEOCEIVER 
STA 31370 1981, and the digital Doppler data, from TRANSIT 
satellite 68, were recorded onto paper tape. Thicty-three 
passes, ovec a period of two andone half weeks, were 
recorded, but analysis of the equipment maintenance 
diagnostics indicated a possible error in the records. 
These data were therefore discarded. 

A second set of forty satellite passes of TRANSIT 
satellite 68 were recorded over a ten day period and were 


determined to be of acceptable quality according to the 
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standards of the DMA Department of Geodesy. This second set 
of data was then submitted to DMA, Department of Geodesy, 
Computation Section, for analysis. Also submitted were 
meteorological data taken during the passes and additional 
information required for documentation of a first order 
TRANSIT position. Satellite 68 Doppler data for the survey 
was then processed in conjunction with the tracked satellite 
poSition by DMA computer facilities. A final geodetic 
position on the WGS-72 ellipsoid was obtained for GEOCEIVER 
STA 31370. Coordinates of the station are listed in the 
Test Procedures section. The computed height relative to 
the ellipsoidal surface is -33.39 meters. Computation of 
the geographic position utilized a minimum of 38 individual 
TRANSIT satellite passes recorded for GEOCEIVER STA 31370. 
This provided for an accuracy of one part in 10° (first 
order accuracy standard). [Ref. 38] 

To orient the positioning network, an initial azimuth 
from GEOCEIVER STA 31370 to station MUSSEL 1932 was computed 
from two sets of observations of POLARIS. All observations 
and computations of the star Sighting utilized procedures 
Outlined. in the National Ocean Survey Photogrammetric 


Instructions Number 4, Revision 1, 4/24/73, the Hydrographic 


89 





Manual, andthe Manual of Geodetic Triangulation, SP-24U7. 
(ref. 39, 40] 

Having the initial azimuth, distances were measured 
between all three control stations using a Tellurometer MRA 
5 microwave distance measuring unit. Twenty observations of 
distance were made for each of the three lines measured: 
GEOCEIVER STA 31370 to MUSSEL, GEOCEIVER STA 31370 to USE 
MON, and MUSSEL to USE MON. The average of the distance for 
each line was then corrected for instrument error and 
meteorological conditions Following the instructions 
contained in the Operators Manual and the Tellurometer 
Manual. [Ref. 41] -Resulting slope- distances were then 
reduced to geodetic distances using the respective 
corrections for slope to horizontal distance, reduction of 
elevation to sea level, and the chord to arc distance. 
These corrections were computed using the WGS 1972 
ellipsoidal parameters and the methods outlined in the 
Hydrographic Manual, the Manual of Geodetic Triangulation, 
and Introduction to Surveying. [Ref. 42] The resulting 
geodetic distances meet third-order, class I specifications 
for traverse accuracy, 1:20,000. In addition, angles to the 


other stations were observed at each control station, giving 





a total of three sets of angles. Each set was comprised of 
four direct and four reverse plate observations with a Wild 
e-2 theodoiite. Each pair of direct and reverse 
observations was corrected for initial pointing errors. 
Observations and computations were made using the methods 
detailed in the Hydrographic Manual, the Manual of Geodetic 
Triangulation, and the Manual of Second and Third-Order 
Triangulation and fTraverse. (Ref. 43) The resulting 
closure meets the specifications for traverse accuracy of 
three seconds per traverse station or a closure of 6 
seconds. 

With the geodetic position of GEOCEIVER STA 31370 known, 
as well as the azimuth from GEOCEIVER STA 31370 to MUSSEL 
and the distance from GEOCEIVER STA 31370 to MUSSEL, the 
gecgraphic position of MUSSEL was computed on the WGS 1972 
ellipsoid using the direct geodetic computation. This 
computation used the Rainsford's method modified by tT. 
Vincenty and including the dHelmert's elliptical terms. 
{Ref. 44} 

Next, using the azimuth from MUSSEL to USE MON, the 
geodetic distance, and the results of the direct geodetic 


computations, the geographic position of USE MON was 





determined. Also, the ARTEMIS geographic position was 
calculated as an offset using an azimuth and distance from 
GROCELVER STA 31370. With these computations, the geodetic 
survey was completed. 

Positioning of the ship, R/V ACANIA, required one 
aZimuth, or line- of-position, from each of the control 
stations MUSSEL, GEOCEIVER STA 31370, and USE MON. These 
lines of position were used to triangulate the vessel 
position. Angles from each of the three theodolite 
locations were observed as described in the Test Procedures. 
Each of the three angular measurements were then corrected 
for mean pointing errors and the individual theodolite 
collimation errors. Each corrected angle was then adjusted, 
by control station, to obtain an aZimuth from south, as 
reguired to triangulate the vessel position. The three 
resulting azimuths and the corresponding control station 
latitude and longitude were then used ina least squares 
three line of position intersection computation. Results of 
the computation were a vessel position in latitude and 
longitude, error ellipse semi-major radius, error ellipse 
semi-minor radius, azimuth from south of semi-major radius, 
and circular standard error. (Ref. 45] Figure 5.1 


illustrates these error ellipse parameters. 


Be 





ViGeeanwm ELLIPSE 


The mean displacement of a point to be expected in the direction of 
the two coordinate axes, X and Y, due to measuring errors, can be ex- 
pressed by the quantities | and iE 


When the north-south axes are rotated, two axes (X and Y) perpendicular 
to one another are obtained. corresponding to the maximum and minimum 
values of © and om . The maximum and minimum values of o, and the 
clockwise rotation atigle (x) of the major axes of the error ellipse are 
included as part of the output of the Geodetic Package and can be used 
as input if desired. 


The circular standard error (CSE) is also included in the output with the 
error ellipse parameters. The value of CES is the radius of a circle 
whose area 1S equal to the area of the error ellipse. That is: 


Coe =| ov Baie 
; max min 





EXAMPLE: Error Ellipse 


Az-Max 144 ,9° 
Max Qi235am 
Min O.21] m 
USE 0.223 m 
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For this study, 262 vessel positions, or fixes, were 
observed and recorded. On examination during processing, 51 
fixes were found to contain errors and data omissions which 
rendered the vessel positioning computation invalid or 
indeterminate. The remaining 211 fixes comprised the data 
set which was then broken down into smaller sets for 
analysis and comparison. These sets were numbered 1 through 
13 and categorized according to vessel track. All sets were 
then processed using several different methods to determine 


the most accurate approach. 


B. THEODOLITE i> Shears SHIP POSITION DETERMINATION 
Before processing the GPS Manpack data, the ARTEMIS data 
were processed to ensure the vessel positions, as obtained 
by theodolite, were correlated with the correct automated 
and hand recorded GPS data. ARTEMIS data consisted of an 
angle, H4USSEL-ARTEMIS OFFSET-vessel, and distance from the 
ARTEMIS OFFSET to the vessel. The angular measurement was 
first converted to an azimuth from south uSing the azimuth 
Of ARTEMIS OFFSET to MUSSEL as the correction and adding the 
measured angle. Vessel position based on ARTEMIS was then 
computed using the direct geodetic program (also used for 


the basic survey}. Comparisons of the resulting geographic 
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positions and the theodolite derived geographic positions, 
were made for corresponding fix numbers using the inverse 
geodetic computation on the WGS-72 ellipsoid. The inverse 
computation utilizes two geographic positions as inputs for 
calculating the forward and back azimuths and the distance 
between the positions. [Ref. 46] For comparison, if the 
two positions differed by more than 20 meters, the ARTEMIS 
positions on either side of the suspect ARTEMIS fix were 
then used to compute the inverse computation ‘to the 
theodolite ship position. In seven fixes, the fix position 
values were changed in the May 15 GPS data to give more 


accurate results. 


C. DIFFERENTIAL GPS SAIP POSITION DETERMINATION 

With the check and adjustment of the GPS data and the 
theodolite determined vessel position complete, the fase 
GPS data processing was the comparison of the uncorrected, 
non-differential ship GPS Manpack receiver data to the 
theodolite computed ship position. Using the inverse 
geodetic computation, the forward azimuth and the distance 
between the theodolite computed ship position and the 
nen-differential GPS Ship position were completed. These 


results were then used in the analysis of the differential 


eS 





processing to examine the stability of the shipboard Manpack 
receiver compared to the GPS Manpack receiver located over 
GEOCEIVER STA 31370. 

There were two methods of analysis of the differential 
GPS. The first of these was the determination of a static 
correction to be applied to all the GPS data. This static 
correction in the differential GPS computations was to 
account for a suspected bias between the two similar, but 
not identical, Manpack receivers. A mean static correction 
was computed for the most stable reception period during the 
static tests on 11 May. With the two receivers collocated 
over the GEOCEIVER STA 31370, a mean correction of -0.954 
seconds of latitude and -0.307 seconds of longitude was 
found between GPS Manpack System 1 and System 2 during times 
07316200 GMT and 07: 25:12 GMT. This static correction was 
then added to the latitude and longitude of GPS Manpack 
System 1 (located over GEOCEIVER STA 31370) data. 

The second method of analysis or processing in the 
differential mode used the latitude and longitude obtained 
directly from the GPS Manpack receivers. In either 
processing method, the computation of the differential 


corrector to the shipboard GPS Manpack receiver data was the 
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same. An explanation of the computations and processing of 
the differential correction for both methods follows. 
Computation of the differentiai corrector for each GPS 
fix position was accomplished by using the geographic 
position obtained with the stationary HManpack receiver, 
located over GEOCEIVER STA 31370, and the geographic 
position of GEOCEIVER STA 31370. These two positions were 
then used in the inverse geodetic computation to obtain a 
differential correction vector (forward azimuth, and 
distance) from the stationary GPS receiver to GEOCEIVER STA 
31370. In figure 5.2, the differential correction vector is 
labelled "A", Using the direct geodetic computation, the 
differential correction vector was applied to the 
corresponding shipboard GPS Manpack receiver position, and 
the differential -GPS position was obtained. These resulting 
positions were computed on a one-to-one basis using the fix 
humbers as identifiers. Vector WB the differential 
corrector vector which is identical to the vector "A" in 
figure 5.2, is applied to the shipboard GPS Manpack receiver 
and results in the differential GPS ship position. 
Comparison of the differential GPS Ship position with 


the theodolite-computed ship position was accomplished using 
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the inverse geodetic computation for both methods of 
differential GPS processing. The resulting forward azimuth 
and distances for both processing methods were used in the 


statistical analysis. 


D. STATISTICAL DATA ANALYSIS 

Statistical analysis was done on data separated into 
sets according to the method of data processing. The static 
data and stability data were graphed as functions of 
latitude and longitude versus tine. User range error, 
furnished by the Space Division, GPS/OL-AQ, Vandenburg, AFB, 
Was plotted as a function of time for each day of the test 
(see Results). Other graphed relationships are the 
gecgraphic positions of the stationary GPS Manpack receiver 
data versus the GEOCEIVER STA 31370 and the theodolite 
semi-major radius as a function of fix number (see Results). 

Other statistical data analyses were represented by 
histograms. Plotted on these histograms are the relative 
frequency on the left margin, variable bar increments scaled 
on the bottom, number of variable observations in each bar 
increment on the top, andthe total sample size on the top 
Margin. Within the histogram the frequency bars are denoted 


by asterisks, the mean by a capital M, the quartiles by 
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periods and the empirical density function by the capital F 


curve. The empirical density function was computed using: 


1 w = 
N(B(N)) i=1 B (N) 


{ Ref. 47} where: F(z) = the empirical density function, z 


the frequency of observation, x = individual observation, 


N = the total number of observations, W = weight function, B 


bandwidth, 
W=1-2 
B(N) = range of all observations / ‘YN 

The statistical values listed below the histogram are 
defined in Table V. Computed histograms show the distances 
between; 1) the theodolite-computed ship position and the 
differential GPS ship position, both with and without the 
static correction, by set; 2) the theodolite-computed ship 
position and the differential GPS ship position, without 
static correction, for five satellite observation periods on 
each of the three ship test days; 3) the theodolite-computed 
Ship position and the non-differential GPS ship position for 


the five satellite-observation periods on each of the three 
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TABLE V. STATISTICAL DEFINITIONS 


The definitions of the basic statistics computed by the program are 
listed below:. 


MEAN: Average of the sample 


MEDIAN: Mid-Value of the sample, if there are an odd number of 
Sample points, or the average of the two middle values 
for an even number of points 


TRIMEAN: 0.25 * (Q, + 20, + Q3), where Q. are the quartiles 


MIDMEAN: The average of all sample values lying between the upper 
and lower quartiles 


MIDRANGE: Average of the maximum and minimum 
GEOMETRIC MEAN: GM ='Nx_-x,-.. . xX 
Nieee2 n 
HARMONIC MEAN: HM = on / 211%, 
VARIANCE: Unbiased estimators for variance and standard deviation 


are used - the square of the standard deviation 
STANDARD DEVIATION: S$ (Fe Ee 
& ater os a ] 


COEFFICIENT OF VARIATION: Standard deviation / |Mean| when the 
mean is less than 10. ** -30, the coeffic- 
lent of variation is set to zero. 


MEAN DEMIATION: The average of the sum of the absolute differences 
between the sample values and the median 


RANGE: Maximum - minimum 

MIDSPREAD: Q3 - 1° also called the interquartile distance 
M3: Third central moment. Unbiased estimator is used 

M4: Fourth central moment. Unbiased estimator is used 
COEFFICIENT OF SKEWNESS: M3 / (Standard Deviation) ** 3 
COEFFICIENT OF KURTOSIS: (M4 / (Standard Deviation) ** 4) - 3 


BETA): Biased estimate of third central moment; can be used in 
testing for normality 


DETAc: Biased estimate of fourth central moment 
MAXIMUM: Largest sample value 
MINIMUM: Smallest sample value 


QUANTILES: The alpha-quantile, x(alpha), is the solution to the 
equation: 
Probability (x .le. x(alpha) ) = alpha 
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ship test days; and 4) the GEOCEIVER STA 31370 position and 


the uncorrected stationary GPS position for the three days 


of ship tests. imeaddition, the azimuths between the 
theodolite-computed static corrections, for each subset, 
were graphed as a series of histcgrams. To illustrate the 


difference between the GEOCEIVER STA 31370 position and the 
stationary GPS position, the aZinuths between the two 
positions are graphed on histograms for each of the three 


test days. 
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VI. ERROR ANALYSIS 








The error sources contributed by each positioning system 
are described. Estimates of the effect of each on the final 


results are determined where possible. 


A. ARTEMIS POSITIONING SYSTEM 

The sources of error in the ARTEMIS ship position can be 
divided into two categories--those sources of error inherent 
in all microwave positioning systems and those unique to the 
ARTEMIS system. All are investigated and the contribution 
of each to the total error is determined. 

1. Microwave System Errors 

Sources of System error include range holes, 

inaccuracies in calibration, non-ideal geometry of the lines 
of position describing the fix, and interference with 


reflections of the signal from nearby objects. 


Range holes are a cause of unreliability in 
Microwave positioning systems. The problem exists when the 
sea surface is calm enough to reflect part of the 


transmitted signal. The reflected signal will arrive at the 
receiving antenna at a later time than the signal which 
travels directly, causing a phase difference between the two 


Signals. A zone or region of destructive interference of 
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the directly transmitted signal and the reflected signal 
ort occur within which the vessel will receive an 
attenuated signal. When the phase difference of the two 
paths is exactly 180°, no signal will be received provided 
the reflected signal is as strong as the direct-path signal. 
Prediction of the range holes, therefore, is useful in 
survey planning so as to avoid this problen. 

The range hole is a function of Shore unit and ship 
unit antenna heights, the wavelength of the system, and the 
distance from the shorebased unit. If the operating area is 
known, antenna heights can be determined and equipment can 
be installed so that no range holes exist in the survey 
area. If fixed antenna heights are known, the region of 
uncertainty can be determined from the equation: [Ref. 48] 


The following equation is used for this computation: 
D, = 2(h1) (h2) /n), 


where A is the wavelength (constant for ARTEMIS at .032 an), 
hi and h2 are antenna heights, and D, is the distance from 
the shore based unit to the nth range hole. For this study, 
the antenna height of the fixed unit above mean sea level at 
the GEOCEIVER STA 31370 was 3.8 meters. The mobile antenna 


Was installed on top of the equipment shelter on the boat 
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deck of the R/V ACANIA at a height of 7.9 meters above the 
the waterline. Two range holes within the survey area were 
determined to be located at ranges of 1876 meters for n= 1 
and 938 meters at n = 2 (fig. 6.1). These two ranges were 
computed for the case of mean lower low tide. The change in 
Width of the range hole due to a tidal fluctuation of 1.5 
meters during the test is shown in figure 6.1. (During the 
Night of 15 May the sea surface was sufficiently disturbed 
by a two-foot swell and a surface wind of 2-4 knots such 
that surface reflection was minimal and no loss of signal 
was observed at the computed ranges). 

Calibration of the ARTEMIS system was accomplished 
in two stages. As part of the study, the ARTEMIS was 
obtained for a one week test of the system from 13 to 17 
April. During this time it was calibrated for azimuth and 
Cange over the baseline between the GEOCEIVER STA 31370 and 
MUSSEL. The azimuth readings required no adjustment but the 
range was reduced by two meters. The equipment was returned 
to Marinav Corporation to fulfill prior commitments and then 
was returned to Monterey for the ship test on 15 May. 
Before the evening operation, the ARTEMIS was used to 


measure the same baseline as opreviously. A plus two meter 
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bias was observed at this time, but no immediate 
compensation was made. This error was assumed constant 
throughout the test and was included in the calculation of 
the ship-determined ARTEMIS position during postprocessing. 

One of the advantages of the ARTEMIS system, unlike 
other microwave positioning systems, is that the geometry is 
always ideal for the lines of position that determine the 
Eo. %. The range-azimuth technique of determining the vessel 
position defines an uncertainty ellipse having axes that are 
always mutually perpendicular. This allows prediction of 
errors as a function of range only, independent of the 
angle. 

The manufacturer's specified error in range 
measurement of the ARTEMIS is + 1.5 meters for all line of 
Sight ranges. One axis of the uncertainty ellipse is 
therefore constant. The other axis, that describes the 
error introduced by the system azimuth uncertainty, is 60 cm 
per 1000 m range. (This is determined from an angular error 
of 0.033 or 2' of arc). Both system error values represent 
an accuracy of twice the standard deviation value, or a 95% 
probability of the true value occurring within the margin 


allowed by these parameters. 
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The positional uncertainty can be expressed in terms 
of the radius of a circle of equivalent probability within 
which the true fix is expected to occur. The method of 
computation of this radius is described in Appendix A. The 
technique allows the conversion of uncertainty ellipses to 
radii of circles of equivalent probability. Values from 
this computation for the radii of circles of 90% probability 
at different ranges used in the survey area for 15 May are 
displayed in figure 6.2 along with the circles of 90% 
probability, the magnitude of the rcadii at the ranges used 
for the study do not exceed 5.4 meters. One goal of this 
study was to measure the difference in positional accuracy 
between ARTEMIS and GPS. A maximum error of 1.0 meters, an 
order of magnitude smaller that the tested system accuracy, 
is desirable. At the present time, however, there are few 
systems that can provide the order of accuracy system for 
the position of a moving vessel, so the 5.4 meter error 
observed at 4000 meters had to suffice. 

Error caused by reflection of the signal from nearby 
objects is difficult to measure. To minimize this problen, 
the offset position for the ARTEMIS unit was planned so as 


to be clear of possible interference from the wooden 
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structure and the eyuipment vehicle. No signal interference 
was observed despite the ship's motion. This source of data 
unreliability was therefore not considered to be 
Significant. 

2s ARTEMIS System Errors 

Errors unique to the ARTEMIS system include a 
pointing error source andthe lack of resolution of the 
azimuth measurement in the dynamic mode of operation. 

The operation of the ARTEMIS system requires a 
manual alignment of the fixed antenna so that it is 
perpendicular toa known object. This is necessary to 
establish an original reference direction upon which all 
angular measurements are based. It is possible that a 
pointing error could be produced by the operator if the 
crosshairs of the telescope used to align the antenna are 
not exactly sighted. This would introduce a systematic 
error into all subsequent angle measurements. To minimize 
this source of error, the alignment of the antenna was 
checked before and after the overiod of operation. The 
Original azimuth displayed onthe antenna unit display 
agreed exactly with the final azimuth vaiue. Any 


contribution to ship position error was considered mininal. 
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The ARTEMIS was operated using the dynamic mode 
option (see Equipment Description--ARTEMIS). The resolution 
of the system in this mode is 1 meter for displayed range 
and 0.1° for azimuth values. This resolution limit 
introduces an uncertainty of 0.5 meters in range and 0.05 
in angle measurements. The computed uncertainty ellipse as 
expressed by radii of circles of equivalent probability, is 
shown in figure 6.2. The system error and all other errors 
mentioned above are assumed to occur within the scope of 
this resolution uncertainty. 

Another source of system error has been proposed in 
a study conducted in June, 1980, by NOAA. (Ref. 49] The 
test compared the ARTEMIS system with two ranges measured by 
two Del Norte Trisponders (a short range microwave 
range/range positioning system) and two lines of direction 
measured by two shore-based theodolites. The results 
described an underway bias of -0.05° of unknown origin in 
the angle measurements. This discrepancy has not been 
resolved. The preliminary test conducted in Monterey in 
April, 1981 compared the ARTEMIS system with three lines of 
direction determined by three theodolites. The purpose of 
the test was to resolve the underway bias problen. Results 


are enclosed in a separate report. [{Ref. 50] 
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B. CONTROL NETWORK 

The entire survey was oriented relative to the TRANSIT 
position on the Coast Guard Pier (GEOCZIVER STA 31370). 
Therefore, errors inherent in the position of this station 
mark affected all results. Assuming that the geographic 
coordinates published by DMA/HTC represented the best 
estimate available of the GEOCEIVER STA 31370, the remaining 
Significant sources of error affecting all results were 
survey errors. The control network for the study was 
determined from the traverse and it was assumed that all 
positions used in the study were oriented relative to each 
other on the WGS-72 ellipsoid. A change in the ellipsoidal 
surface within the area encompassed by the traverse was 
assumed to he insignificant. 

Survey errors stem from the methods used for obtaining 
the initial azimuth for orientation of the survey, the type 
of «traverse chosen, and the geometry of the traversed 
triangle. 

To obtain a starting azimuth for the traverse, two sets 
of observations on Polaris were recorded at the GEOCEIVER 
sea 31370. A striding level attachment to the theodolite 


Was not used. The instrument level bubble was observed 
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repeatedly during the observations to ensure constant 
instrument alignment with the vertical. Failure to keep the 
instrument level could have introduced a +t 30" uncertainty 
into the starting azimuth of the survey. [ Ref. 51) To 
consider the significance of such an error, all stations 
were recomputed using a deviation of + 15" from the computed 
azimuth. The resulting shift in the positions of MUSSEL and 


USE MON are listed: 


USE MON 0.133 m distance at 230°36'47.16" azimuth 


MUSSEL 0.132 m distance at 230 33 4.08" azimuth 


The change is considered insignificant relative to the other 
errors described above. 

The survey was a closed loop traverse using GEOCEIVER 
STA 31370 as the origin and closing position. The weakness 
in this method involves the measured distances. A method of 
checking the observed angles in the triangle exists via the 
angular closure if a significant deviation in one of the 
angles can be detected. The distance mneasurement, however, 
Can produce an undetected bias in the survey if the distance 
measuring instrument contains a systematic error. This 


error would shift all the surveyed positions. Corneet 
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measurement of distance by the MRA V Teliurometer, 
therefore, was important. The equipment was calibrated over 
a known Ssetacs and the errors were within the centimeter 
range. This error was therefore negligible. 

The geometry of the triangle was not optinun. Ideally, 
the triangle would have been equilateral +o minimize the 
source of observing error. The choice of traversed points 
was based on the location of preexisting station marks and 
the desired geometry of the theodolite fixes to the ship. 
The closures computed for the traverse met third order 
specifications for accuracy, and the survey was not adjusted 
to remove the closure error. [Ref. 52, 53] 

The survey errors were considered to be within the 
uncertainty cf the TRANSIT station location and were 


therefore not considered significant for the study. 


C. THEODOLITE NETWORK 
The errors generated by the theodolite network are 
divided into equipment errors, Operator related errors and 
System errors, of which the latter is the most significant. 
The equipment used for the theodolite network were three 
1-second theodolites which have a measurement resolution of 


Om". All theodolites were in good conddstion. 
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Instrumental errors of a theodolite are listed in Table VI. 
(Ref. 54] 

The major source of equipment error was introduced by 
not changing face, the process of observing angles a second 
time with the telescope reversed. Due to the nature of the 
test coupled with 11-minute fix intervals, it was not 
possible to observe all the fixes in both the direct and 
reversed modes. AS an example of the magnitude of this 
error, the horizontal collimation errors were determined for 


each theodolite prior to the study. The results are listed 


below: 
Theodolite Collimation gError 
# 14405 ou 
# 14482 cea 
# 14452 i) 7 


All the measured angles were adjusted for this error as 
described in the processing section. The largest influence 
on the remaining instrumental errors was the non-verticality 
of the main axis. If the non-verticality were large (for 
example, one division off-level in the lower spirit level 
bubble), the angle measured could be in error by 7". 


(Ref. 55] To minimize this error, each theodolite was 
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TABLE VI. Theodolite Instrument Errors 


ERROR ELIMINATED BY: 
Eccentricity, misplaced indices Reading both sides of circle 
and changing face 
Graduation of circles Changing zeros 
Micrometer run and vernier error Selection of suitable zero 
setting 
Non-verticality of main axis Error is not cancelled by any 
observing procedure 
Bad setting of the upper spirit Correct setting or applica- 
level GOmeOr COmrection 
Verticality of cross wires Changing face and using same 


part of wire 
Horizontal collimation Changing face 


Vertical collimation and index Changing face 
setting error 


116 


ays 





relevelled if necessary between sets of observations. The 
magnitude of this and the collimation errors were considered 
to be insignificant relative to the system errors discussed 
later in this section. 

The operator-related errors are a function of the 
individual performance of the observers. The experience and 
personal comfort of the observer, for example, are two najor 
factors which can influence the observed angles. An 
observer's performance and resultant impact on the quality 
of the results depend on the judgment and reaction time of 
the observer and are influenced by variables such as fatigue 
and age. The subjective nature of the problem makes an 
estimate of the error difficult, therefore, for this study 
the error was considered randon. 

The major source of error in the theodolite positioning 
system dealt with the timing of the fix. A delay existed 
between the time the ‘mark’ was called over the radio and 
the actual time the angle was observed. The timing problem 
Was most critical when the ship was close to one of the 
three observer locations because the rate of change of the 
bearing of the ship was greatest at this time. Hence a 


Small delay generated a large amount of error in the angle 
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recorded. To indicate the effect this error source may have 
had on the position of the ship, two fixes were chosen 
arbitrarily--one at close range to each of MUSSEL and 
GEOCEIVER STA 31370. A one minute value was added to the 
angle which was most affected by the timing error and the 
fixes were recomputed. The same process was repeated for 
the three fixes with the one minute value subtracted fron 
the critical angle in each fix. The results follow: 
Number Original Pix + 1 Minute - 1 Minute 
210 36°37" 20.62282" 36°37°20.61468" 36°37'20.63081" 
121 53°41.97184" 121 53'41.96933" 121 53°41.97426" 
259 36° 36°45,.88048" 36° 36945.87574" 36°36'45., 89082" 
127 52°54.05062" 121. 52*54.04635" 121.52°54.05459" 
No Significant error was generated from this source. 

Another source of system error was associated with the 
geometry of the fix. Ideally, the three theodolite angles 
would describe two equilateral triangles to the ship; 
however, within the operating area, the geometry varied. As 
an indication of the size of this error, representative 
fixes were chosen and the circular map accuracy standards 
(CMAS) of each of the combinations was determined. Figure 


6.3 illustrates the approximate error caused by variations 
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in the geometry. Within the area outlined in the figure the 
gecmetry creates a CMAS of not more than 10.4 meters. This 
is therefore the source of the largest probable error in the 
theodolite method of positioning. Nevertheless, this method 
in combination with the ARTEMIS system provided the best 


estimate of the ship location. 


D. TRANSIT POSITIONING SYSTEM 

The method used to determine the TRANSIT position 
(GEOCEIVER STA 31370) is the most accurate satellite 
positioning technique presently available. Error sources 
are reduced both during data collection and during 
postprocessing of the recorded values. 

DMA places requirements on the collection of data that 
Will minimize error. For this study, only satellite passes 
whose maximum altitude ranged between 10° and 80° were 
recorded. Meteorological data (temperature, pressure, and 
humidity) were collected for each pass, anda determination 
of the receiver clock drift was computed after each pass. 
The amount of clock drift asa function of time gave a 
real-time indication of receiver performance. Equipment 
checks were conducted routinely, and the results were 


included with the data package submitted to DMA. 
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An error budget of the TRANSIT system for point 
positioning is listed in Table VII. [Ref. 56] The precise 
ephemeris method of postprocessing accounts for all of these 
errors: 1) for each pass the tropospheric correction is 
computed from the meteorological data; 2) the ionospheric 
correction is determined by observing the shift in frequency 
between the two transmitted frequencies; 3) the receiver 
clock error is known; 4) the satellite orbits are observed 
by the tracking stations so that precise satellite positions 
are known; 5) errors due to uncertainties in the 
geopotential model and surface drag forces are accounted for 
during computation; and 6) the receiver altitude relative to 
the ellipsoidal surface is computed approximately from the 
given MSL elevation. 

The accuracy of the computed TRANSIT position is stated 
as 1.5 meters per axis within a 90% confidence interval. 
This is assumed to provide the best estimate possible of the 
GEOCEIVER STA 31370 relative to the wWGS-72 coordinate 


system. 


Ee GPS ERROR SOURCES 
Global Positioning System error sources are separated 


into two main categories, the space vehicle system errors 
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TABLE YII.TRANSIT System Errors 


SOURCE ERROR (m) 


Uncorrected propagation effects 
(ionospheric and tropospheric effects) 1-5 


Instrumentation and measurement noise 
(local and satellite oscillator phase 
jitter, navigator's clock error) 3-6 


Uncertainties in the geopotential model 
used in generating the orbit 10-20 


Incorrectly modelled surface forces 
(drag and radiation pressure acting on 
the satellites during the extrapolation 


interval ) 10-25 
Uncertainties in navigator's altitude 
(results in a bias in longitude) 10 


Ephemeris rounding error 5 





and the GPS user receiver errors. There are seven possible 


sources of error in the space vehicle system--errors 


associated with satellite equipment, satellite clock, 
ephemeris, atmospheric delays, Signal multipath effect, 
satellite update, and satellite propulsion firing. Table 


III lists these error sources and the approximate user 
equivalent range error for each. [Ref. 57] 

The satellite equipment error, also termed the group 
delay, is caused by differences and uncertainties in 
processing and transmission of the message through the 
individual satellite's circuitry. Bach satellite has the 
group delay calibrated prior to launch and this Known group 
delay is included in the satellite time corrector within the 
navigation message. However, unknown delays in signal 
processing and circuitry not accounted for by the corrector 
contributes an estimated one meter to the error budget. 

Satellite clock error is caused by variations in each 
satellite clock. Satellite clock time may depart by up to 
976 microseconds from the correct GPS time. Corrections for 
the predicted clock variations are included in the 
Navigation message. Residual clock variations contribute to 


the error budget and are included in ephemeris errors since 
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the clock errors are small and have the same effect as 
ephemeris errors. 

Ephemeris errors are the result of satellite clock 
variations, drift by the monitor station clocks, monitor 
station signal processing delays, lack of precise monitor 
Station geodetic positions, solar wind and pressure 
parameter variations, and earth gravity model error. These 
sources all contribute to the ephemeris error budget. The 
ephemeris error budget is monitored over long time periods 
and translated by the Master Control Station into an 
User-Equivalent Range Error (URE). When the URE exceeds the 
prescribed value of 4 meters the satellite ephemeris data is 
uploaded to reduce the URE. { Ref. 58] Most of the URE is 
included into the user clock corrector obtained in the GPS 
position computation. However, not all of the ephemeris 
error is removed during the position computation, and the 
residual error has been estimated to be approximately 1.5 
meters. See Table III. 

Atmospheric delays in Signal propagation are the result 
of two effects, ionospheric delay and tropospheric delay. 
Of these, the ionospheric delay, which depends on the 


frequency of the transmitted signal, is greater. With the 
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two frequencies, L1 and L2, transmitted by each satellite to 
the user receiver, the ionospheric propagation delay can be 
modelled. Using this medel, the majority of the propagation 
delay is accounted for, but a residual delay results in an 
error. The other part of the atmospheric delay error 
depends on the tropospheric propagation delay. Propagation 
delay in the troposphere is a function of Signal path length 
through the layer, the receiver air temperature, and air 
pressure and humidity. Since the tropospheric delay is not 
a function of frequency and is highly localized, the 
tropospheric delay corrector has not been incorporated into 
the atmospheric correction computation and is therefore an 
additional source of error. These two error sources combine 
to yield an atmospheric delay error estimated t2 have a 
Magnitude of between 2.4 to 5.2 meters. 

Ionospheric scintillation of the GPS signal, in addition 
to the ionospheric propagation delay, Cause errors in 
positioning. Scintillation has the effect of random 
interference with the propagated signal as the signal is 
transmitted through the upper ionosphere. (Ref. 59] 
Interference caused by scintillation combines with the GPS 


user receiver processing noise to cause a random tranging 
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error and, if the interference is severe, interruption of 
the pseudo-random number binary code causes loss of signal. 
Errors and signal losses caused by scintillation cannot be 
predicted or estimated. 

Signal multipath errors result from the reception of a 
satellite's signal transmitted via different propagation 
paths. Received signals thus may have been distorted or 
reduced in strength to a sufficient degree so as to affect 
the range determination from the satellite. Error due to 
multipathing is highly localized and cannot be modelled. 
Estimates of the error magnitude are between 1.2 and 2.7 
meters. 

The last two space vehicle system error sources are 
attributable to the Master Control Station (MCS) satellite 
accuracy and maintenance functions. Usually the MCS uploads 
each satellite at a given time each day and the ephemeris 
data is changed(Table VIII). Should the user not be aware of 


occurrence of upload and if the receiver therefore does not 


obtain the new ephemeris data, Tanging to the satellite in 
question will contain a possibly very large error. In 
addition, if the més, through the Remote Stations, 


determines that the four skewed reaction wheels cannot 
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Maintain satellite antenna alignment With earth, the 
hydrazine propulsion system is activated and the satellite 
1s rotated. Once the propulsion system has used, the 
satellite orbit cannot stabilize immediately, and the 
ephemeris data will be inaccurate. A period of up to 
several hours may be required for orbit stabilization and 
for the new ephemeris data to be computed and uploaded. 
During that period, the ranging computations to the 
satellite will have large errors. The magnitude of these 
twe error sources cannot be eStimated but is assumed to be 
several orders of magnitude larger than the previously 
stated errors. 

The second error category is that related to the GPS 
user receivers. These sources of error associated with the 
user GPS receiver, the Manpack, are receiver signal noise, 
receiver vessel velocity, Signal resolution, method of 
Signal processing, ephemeris update capability, and, for the 
differential GPS computations, system accuracy comparison. 

Received signal noise errors result from interference 
introduced into the receiver GPS signals through processing 
and internal receiver circuitry. Additionally, quality of 


Signal resolution is dependent upon receiver type and can be 
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a source of error. The combined errors from these two 
sources are estimated to be 1.5 meters. [Ref. 60] 
Estimates of these error sources for the MVUE have not been 
made but are estimated to be of the same magnitude. 

Velocity of the receiver vessel results in GPS errors of 
varying magnitude. The higher the velocity, the larger the 
error. This error is due to the time lag required to obtain 
and process each NAVSTAR satellite range with the MVUE 
Single channel receiver, to compute the average vessel 
velocity and to use the average velocity for positional 
computations. The Manpack requires a six second time 
interval to obtain and process the pseudo-range from each 
satellite used. Thus the fix position is an extrapolation 
of the four satellite ranges ‘+0 obtain a convergent fix 
position. At one knot, the vessel travels 0.51 meters in 
one second. During the full twenty-four seconds required to 
obtain four satellite ranges, the vessel would have 
travelled 12.24 meters. With the vessel travelling at four 
knets, the distance travelled during the total time interval 
required for a GPS £ix position would be 48.96 meters. 
Processing the pseudo-ranges from four satellites requires 


the solution including the average vessel velocity during 
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the 24-second ranging period. Use of this average neglects 
the relative velocity of the receiver antenna which is 
mounted on the mast. During this test, the GPS antenna was 
mounted 9.06 meters above the waterline of the R/V ACANTA. 
Assuming five degrees of pitch and roll anda vessel roll 
period cf 15 seconds, relative average antenna velocity was 
0.16 meters per second. The antenna velocity error 
translates into a random error of 0.16 meters. The other 
velocity errors have not been estimated but are assumed to 
be less than one meter for the low vessel velocities used in 
the test. 

One error source mentioned previously requires further 
clarifications: error due to the relative position at the 
time of the fix. The Manpack requires up to twenty-four 
seconds to obtain and process the satellite ranges. The 
position of the fix and the associated time of fix are then 
recorded. However, the time of the fix may not have been 
the same time the computational results were printed. Thus 
the output position may not have been correlated with the 
correct theodolite computed ship position since the 
theodolite lines of position were observed when the GPS 


position was printed. Errors from this bias source may 
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cange from 12.24 meters at one knot of vessel velocity to 
48.96 meters at four knots. This error varied with the 
precessing time required to obtain a convergent position. 
An estimate pf Miktis error magnitude was not made because of 
its wide variability. 

The type of receiver and the method of signal processing 
also contributed to the error. Multichannel receivers 
designed for Phase I testing have a lower magnitude error 
for satellite signal reception and processing than do the 
Single channel receivers, such as the Manpack. Ine *ehe 
multichannel receivers, four or more satellite pseudo-range 
data are taken simultaneously. The single channel receiver 
seguentially receives each satellite and computes the 
pseudo-range data. Simultaneous data reception reduces the 
user velocity error and therefore has less position aliasing 
than that found in the single channel receivers. 

The method of signal processing may induce positional 
errors through the use of non-optimal weighting functions in 
the position computation. Should the ephemeris or ranging 
data change significantly, as happens during uploading, 
Signal processing using Kalman filters as used in the 


Yanpack do not change the computed position rapidly enough 


13] 





to reflect the large changes in ranges and result in a slow 
change drift of the GPS position. EStimated tine 
requirements for user receivers to stabilize after an upload 
are up to two hours. [Ref. 61] 

The capability of the user GPS receiver to receive 
ephemeris data immediately after upload reduces the 
ephemeris errors caused by using an outdated ephemeris. The 
ephemeris data is repeated every thirty seconds and its 
length requires sections of two message subframes. 
(Ref. 62, 63). The use of current ephemeris data in 
the computations reduces the error for the ephemeris to the 
value shown in table III. The Manpack has the capability to 
accept new ephemeris data immediately, but during the test 
large errors were introduced into the computation on 
satellite ephemeris upload because this selected optional 
function either did not operate properly or was neglected by 
the operators. 

For the differential GPS user, an additional error 
source is introduced into the error budget: the assumption 
that both GPS receivers will observe, obtain the 
pseudo-ranges, and compute the position solution in exactly 


the same manner. If <+¢his were true then both receivers 
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Should be able to compute the identical position if placed 
at the same position at the same time. With the Manpacks, 
the assumption that the units were identical and the 
position solutions were identical was not made because prior 
testing established that the receivers differed in position 
computations by as much as 3.08 meters. [Ref. 64] 

For the static testing on 11 May, one error source not 
previously discussed was the vertical separation of the two 
Manpack antennas, colocated over the GEOCEIVER STA 31370 
Geutsol station, and the possible interference of the 
uppermost antenna, System 1, with the GPS signals received 
by the lower System 2 antenna. Separation of the two 
antennas was 2.1 meters. Calculations using the speed of 
propagation in avacuum and the two GPS frequencies produce 
wave lengths of 0.19 meters and Q.24 meters. Since 
separation of the antennas was much greater than either 
wavelength theoretically there should have been no 
positional errors due to antenna interaction or 
interference. Differing elevations of the antennas above 
the ellipsoidal surface also may have induced errors in the 
position computation. Estimated errors for the elevation 


differences are less than 0.1 meters. Total errors for the 
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antenna vertical elevation separation and antenna 
interference or interaction have not been estinated. 

Thus, the error budget for the Global Positioning System 
Suggests a combined error of up to 57.35 meters. Using the 
differential technique, the bias errors can be reduced in 
Magnitude or eliminated. These bias errors include 
ionospheric propagation delays, tropospheric propagation 
delays, variations from predicted ephemerides, and satellite 
clock perturbations. Teasley, Hoover, and Johnson (1980), 
evaluated the error budget for a more advanced GPS receiver 
and found the predicted filtered solution error was reduced 
from approximately 11 meters to approximately 3 meters with 
differential GPS processing. fRef. 65] Error budgets for 
the Manpack and the differential correction vector 
processing method used in this test did not attain the four 
meter accuracy but were not expected to because the 


receivers used were less accurate, and of different types. 
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Vit. BESTS 


Results of the GPS tests have been separated into four 
sections: stationary GPS evaluation, evaluation of non- 
differential GPS, differential GPS evaluation with static 
correction, and evaluation of differential GPS without 
static correction. Tables IV and VIII list the optinun 
satellite observation periods and the NAVSTAR upload times 


for each day of the testing period. 


A. STATIONARY GPS EVALUATION 

Stationary GPS receiver evaluation was separated into 
the stability tests, the static tests and the tests of the 
stationary GPS receiver used during differential testing. 
Figures 7.1, 7.2, and 7.3 illustrate the results for the 
stability test conducted on 9 and 10 May. The receiver was 
less stable than expected and depended to a great extent on 
the NAVSTAR uploading times and the URE of NAVSTAR 1 (fig. 
7.4 and 7.5). Correlation between times for the NAVSTAR 
upload and changes in the stability can be observed at 0620 
Oa 9 Bay (fig. 7.1), 0655 on 9 May (fig. 7.2), 0815 on 10 
May (fig. 7.8), amd OS4S on 10 May (fig. 7.3). The 
relationship between stability and URE is illustrated by 


comparison of the JRE for NAVSTAR 1 at 0640 on 9 May (fig. 
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Peay with the stability plot on figure 7.1 for the same time 
period. On figure 7.2 the stability plot variation at 0745 
on 10 May remains unexplained but may have been due toa 
receiver malfunction. 

Results for tests on 11 May (fig. 7.6, 7.7, and 7.8), 
exhibited less stability than expected. Also, system 2 was 
less stable than system 1. The largest difference in 
reading between the two receivers occurred at 0750 (fig. 
wT). This difference may have been caused by the location 
of the antenna of system 2 under the antenna of system 1. 
At approximately 0800, the satellite constellation elevation 
angles were the largest, between 56 to 80° above the 
horizon (fig. 7.9), and antenna interference may have led to 
the large differences. Neither the satellite constellation 
aZimuth angles (fig. 7.10) nor the user range errors Cand’. 
7.11) appeared to have any effect on the stability. 

Comparison of the stationary GPS receiver positions with 
the position obtained with the geoceiver is depicted in 
figure 7.12 for the stability and static tests. The 
comparison indicates that there ls a nearly constant offset 
of the two positions with the GPS position almost due south 


of the geoceiver position. Further comparisons between the 
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stationary GPS receiver position and the geoceiver pcesition 
reinforce this conclusion. In figures 7.13, 7.14, and 7.15, 
the mean distances between the two positions are 127, 173 
and 143 meters, respectively. The higher means for the data 
shown in figures 7.14 and 7.15 are due to the URE (fig. 7.16 
and 7.17) for those dates and the NAVSTAR upload frequency. 
If the data from 13 May is deleted because of the high URE 
for that date (fig. 7.16), and the median distance used for 
15 May is taken to be a representative value (130 meters on 
Fig. 7.15), The difference is approximately 128 meters 
Between the two positions. The azimuth was also nearly 
constant with the observed median values of 181, 207° and 
178° from south as shown in figures 7.18, 7.19, and 7.20 
respectively. Neglecting the data on 13 May because of URE 
values (fig. 7.16), the remaining two values are nearly 
identical and have a mean of 179.5 from south. Thus the 
data indicates a nearly constant offset between «he position 
Obtained with the geoceiver and the position obtained with 
GPS. This difference is about 130 meters and has an almost 


north-south orientation. 
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Figure 7.14. 
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Figure 7.20. 
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Be. EVALUATION OF NONDIPFERENTIAL GPS 

Evaluation of the non-differential GPS positions was 
done on a set by set basis. Figures 7.21 through 7.31 depict 
the distance in meters between the theodolite-computed ship 
position and the shipboard GPS receiver position. Sets 1, 2 
and 3 have the smallest distance difference, as shown in 
figures 7.21, 7.22, and 7.23, with mean values between 118 
and 133 meters and standard deviations ranging from 7.45 to 
16 meters. These smaller mean distances are a result of the 
low user range errors for 12 May as shown on figure 7.32. 
The highest mean distances occurred during sets 4, 9, 10, 
12, and 13 when the JURE values were high and the NAVSTAR 
uploads occurred (fig. 7.16 and 7.17). The remaining 
distances vary between 144 and 278 meters. 

The results during the best satellite observation 
periods are depicted in figures 7.33, 7.34, and 7.35, where 
the non-differential GPS ship position was compared to the 
theodolite-computed ship position, the distance between the 
two being plotted as a histogram. Mean distances on figures 
7.33 and 7.35 are representative of the five satellites 
observed and the improved accuracy of the system, with mean 


distances of 134 and 149 meters and standard deviations of 
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Figure 7.2]. 
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Figure 7.33. 
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10 and 52 meters. The large standard deviation for 15 May 
is apparent in the bimodal empirical frequency curve (fid. 
7.35) and is caused by the NAVSTAR upload (fig. 7.17). In 
this case, the GPS receiver did not obtain the new ephemeris 
until halfway through the time period of the observaticns. 
The data for 13 May (fig. 7.34%) reflects frequent NAVSTAR 1 
uploads and the resulting high URE values (fig. 7.16) for 


that date. 


C. DIFFERENTIAL GPS EVALUATION WITH STATIC CORRECTION 

The results of the evaluation of the data of 12 May 
(sets 1, 2, and 3), computed using the differential GPS 
positions and the static corrections, are shown in figures 
7.36, 7.37, and 7.38. Mean difference between positions for 
the theodolite ship position and the GPS position are 29.1, 
27.4, and 64.8 meters, with standard deviations of 7.7, 7.5, 
and 22.7 meters. These distances are larger than those 
determined for the theodolite ship position and the 
differential GPS position (fig. 7.39, 7.40 and 7.41). This 
higher mean distance indicates that the static correction, 
as applied, results in a larger error in the differential 
GPS position. Therefore, further evaluation of differential 


GPS with static corrections were not conducted. 
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Figure 7.37. 
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Figure 7.39. 
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Figure 7.40. 
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D. EVALUATION OF DIFFERENTIAL GPS WITHOOT SPALIC 
CORRECTIONS 

Evaluation of differential GPS wacthowt the static 
corrections were computed by sets for the five satellite 
constellations for each day of testing. Set 1@@eig. 7239) 
Shows a mean distance of 10.1 meters between the two 
positions anda standard deviation of 5.0 meters. The 
median azimuth between the two was 257° from south (fig. 
7.42), which also was the approximate course of the vessel 
during the test. These values indicate that the error is 
caused by the bias, or alias, in GPS position due to the 
Manpack 24-second fixing interval and positioning algorithm. 
For this set, the GPS position was obtained after the 
theodolite position. 

During set 2 the mean distance difference between the 
two positions was 11.7 meters, and the standard deviation 
was 6.8 meters (fig. 7.40); the median azimuth was 77 
Berg. 7.43). These values were slightly greater than those 
of set 1 due to the increased URE for NAVSTAR 1 during the 
set (see 0700 on fig. 7.32) Additionally, the azinuth was 
nearly 180° from that computed in set 1, with the vessel on 


the same track as set 1. Onlike the error obtained from the 
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Figure 7.42. 
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Figure 7.43. 
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Manpack bias onset 1, the error in set 2 was due toa 
Similar time delay bias because the theodolite position was 
obtained after the GPS position. 

Set 3 values on figures 7.41 and 7.44 are similar to 
those of set 2, except the mean distance values were 
greater--44.7 meters for set 3 versus 11.7 meters for set 2. 
Mean distance increases for set 3 were due tc the upload of 
NAVSTAR 1 (fig. 7.32) and the availability of only a four 
satellite constellation. Set 2 had a five satellite 
constellation available. 

For set 4, figures 7.45 and 7.46 show the distances and 
azimuths for the set. These histograms fail to illustrate 
the constantly increasing distances between the two 
positions. These increasing distances result from the 
uploads of NAVSTAR's 5 and 6 and the large URE for NAVSTAR 1 
feag. 7.16). Onboard the vessel the new ephemeris data 
uploaded for NAVSTAR‘*s 5 and 6 was not obtained until late 
in the set. Also near the end of the set, NAVSTAR 1 was 
uploaded with the result that the last several distances 
were very large, as shown in figure 7.16. 

En )set 95, a very high URE value for NAVSTAR 1 


contributed to the extremely large position difference 
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Figure 7.44, 
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Figure 7.45, 
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Figure 7.46. 
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between that determined by the theodolites and that for the 
differential GPS position (fig. 7.47). The azimuth remained 
nearly constant between the two positions with a median 
aZimuth of 43. (fig. 7.48). These values for set 5 show 
that the large URE for NAVSTAR 1 was the most important 
contributing error source (fig. 7.16). The large URE may 
have been caused by activation of the hydrazine 
stabilization jets to dampen momentum on NAVSTAR 1. 

For set 6, the mean distance difference between the two 
positions was 60.0 meters and a median azimuth of 93. (Eg). 
7.49 and 7.50). These values indicate that the upload of 
NAVSTAR 1 near the beginning of the set (0715 on fig. 7.16) 
was received and utilized in the Manpack position algoritha, 
but the higher distance between the positions also indicates 
a large URE (also shown on fig. 7.16). 

Data from set 7 were not evaluated because the 
Stationary Manpack receiver did not record the GPS position 
for each fix due to radio failure. The mean distance 
between the two positions for set 8 was 121 meters with a 
mean azimuth of 68 Giaige 7eo\  gmpd 7/252). Both the 
distances and azimuths show that the URE for NAVSTAR 1 at 


M590 (fig. 7.17) Was increasing in magnitude and thus 
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Figure 7.47. 
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Figure 7.49. 
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Figure 7.50. 
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Figure 7.51. 
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resulting in a greater positional difference. Another 
factor was the bias introduced by the Manpack 24-second 
reception period. This is evident in the azimuth for the 
theodolite ship position obtained from the GPS shipbcard 
receiver. 

The data for sets 9 and 10 were combined for evaluation 
and the results were a mean distance between positions of 
241 meters and an azimuth of 267 (fig, 7.53 and 7.3%). 
Values of that magnitude show that the URE for NAVSTAR 1 was 
the error source. Data were sampled from 0600 to 0640 (fig. 
7.17). Thus the correlation between the high NAVSTAR 1 URE 
and the high mean distance difference between positicns. 
Again, receiver bias also introduced an error into the 
evaluation. 

Analysis of set 11 shows a bimodal distribution in both 
the two position differences and azimuths (fig. 7.55 and 
7 a6). The peaks on the distance histogram are 
approximately 133 and 160 meters; the peaks on the azimuth 
histogram, approximately 220 °and 236°. Both peaks are 
caused by the upload of NAVSTAR 3 halfway through the set 
and the increasing magnitude of the NAVSTAR 1 URE just after 
the upload as shown by figure 7.17. The total positional 


difference for set 11 was approximately 154 meters. 


eS 





Figure 7.53. 


°0 


CUSSESIIOL °E WNWIXYW CO*IRTSBIE°ZS NVIW WHYH 
Eaeaea.- s JTIINYVNO 06° LOR Eee et zvi38 EL ct OV3INdSAIW Slee adie <s NvjwW WI39 
043120402 ° ye rau JVJINVNOD SL° $0439SCCBO0°L Tv¥i39 2043260940 °2 JONTY s31t9elg9'2 JONVHIITH 
Bee oc Ges (NVIQIW) JIILNVNDS Eee eee t sd See ade TO*93ECZIOZL°E A30 NV3W Se aroun: 6 NY IKOTW 
O+38E0C40° (JONIH) JIIANVANOD Se° Ba C%6STI°9 SINMINS 10-36S0E4I°2 ¥VA 4309 943G64€E105 °2 Nv IJbTUL 

C0 *#31LS959°T JVIANVNO O° LO43C9RL09°2Z oW TO*3SE8G RIS AIO Cis CO*ZTITAL 9° 2 Nv JI3n 

20431156653 ° 1 WOWINIW %0#39215S2S°8 tW €0+3182689°2 JONVIUVA 2O*3/2R6T45°? hy 2a 

NOTLNGIYLsSIa SIN3JnOW THINS) YIHOTH OQvauds ADNAONIL WHI) 
SYOLIINNOD LNOHLIM NCTLISOd $dD IWIINSHIISSIT 3H1 ONY 
NOLPTLISOd A@IHS 2b NISMISG SYILIW NI SSINVISIO OT ONY 6 S145 
€ °9S€ "THE "92E *21¢ "162 "282 °892 "€S$2 "6€2 "422 *6 °G5l “cel °991 
I { f | i j i j i j i i 

OOF nt na 8m Ot em Fm an gn tm $$ 68 eH mH 8 HH HF HHH + OOH HK HOH HH EHF EH HOE COHE HE 4H OF H—-4— 9-4 ——-F OHHH KH -—FeHEtE ES 

oe¢ ee9 eee ° #200 wee #00 “he ee e466 448 . eee eee O06 

eee ee6 eee -" 99¢ oe 060 “He @e0e ee Fe8 ce eee oe G08 

$e¢ ¢e% eee " @re e600 ee “he ene e090 e089 a eee eee eee 

ees 3$434344434444354544 e¢ eee °" gee wee eee “be ee e680 208 : eee eee eee 

eee J444 44 e+e eo? " 900 e600 vee “he 6o8 686 O08 : eee eee e¢ 

3444 3j eee eee - ee 660 oh “Ke #60 #60 e080 , eee ee eee 

en¢ 4-4 eee eee * @66 060 eee “Ne 0 vee e8 : eee eee eee 

oe6 J eee eee ° @86 266 60 “he e860 e080 ve8 : oe% oe etel20° 

ees J ee aoe " eee 206 e9e "Ke 266 200 O28 : een oon oF 

ete : sSjieoe ee¢ - @¢e 66 #260 “He Coe 206 oF . eee eee 460 

eee 4ee eee * @e¢ 00 e680 “Ne ser 66 £68 : eee eee eo@ 

eo jie see ° e¢¢ eee ate “he 0 eee 268 : eee eee ¢98 

eee eed ee6 " ee tee see “he 606 e006 ve% - eee eee 968 

+o0 eee eee " 600 60 ¢¢¢ “Ke @6e te G28 . eee eee bee 

ee% eoeds ee6 * eee 200 ee0 “We eee eee e080 , ees eee e¢8 

eee eee 3 eee " 00 tee eee “he Gee tee ef8 i eon eee ee 

eve eee 3 ¢e8 * 60 eet e880 “Ke Fe O08 e648 : eee eee e2¢150° 

+e ee 3 eee * 666 eee ee “He 666 eee oer : eee Cee G06 

eee ee6 4 eee " 60 @f6 00 “We wee ote een a eee eee eee 

ee8 ee ¢ 4 eee * 860 @2¢ ee “he eee eee ee8 . e6 eee ere 

e+e ¢¢0 4 eee " @06 eee e840 “he e086 60 e680 ° 6 eee ee4 

169 eee oe " @96 @0¢ e668 “ht O06 O06 e086 * eee J4ttjedgisiea 

eee : eee d eee ° 660 006 06 “We eee ee 4494 c eed4 ese eee 

ose eee 4 oe8 “ @te e66 eee “he 60 608 #08 : 4f4iie eee eee 

3 e689 : een “WH 26% "43 ¢e2190° 
ee¢ . eee °H eee © 4 ee@ 
4 e080 . eee “WH 68 jj eee 
eee 7 eee “WN @ee 4° eee 
4906 a eee “WH e806 jj ° eee 
see : eee “WH ee . eee 
soe = eee “HN e086 4 : eee 
ed : eee “WN eee 3 : ee 

ee > eee “WH @ee J 3 eee180° 
eee5 °* oe “NWN oe 3 . ee 
ee¢ : eee “WH eee : eve 
eee Jj ° eee “WH e6 J . ee6 
eee 4° eee “WH oe 3 ‘ eee 
ees 3 nes “WH e¢@ 4 z ene 
eee °4 eet “WN 080 J j eee 
eee = eee “WH tes _ eee 
eee °* J eee “W eee J i eee 

: eee ° 3 ooo “WH eee Jf : eee sO? 
ee6 . 3 eee “WN 48 2 eo 
e¢ e eee “WN eee = e 
oes 33 eee “WH o04 . eee 
ee9 a Jj eve “WH ede : ene 
eee = dooe "W 4t4e0 : eon 
eee - Sie “Jd eee i eee 
see : eeddijWN e0¢@ ° eee 


1 0 6—O 0 0 0 oO 9 1 0.600 0 6 oe | t z 1 Z T t o oO 0 1 2 1 2 
LT = 3218S JVWVS SITININS3I wd 


196 





Figure 7.54. 
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Evaluation of the combined sets 12 and 13 (fig. 7.57 and 
7.58) result in a mean distance between the two positions of 
172 meters and a mean azimuth of 262°. The mean distance 
was large because of the large magnitude of the NAVSTAR 1 
OMe (fig. 7.17). A contributing source of error was the 
bias introduced by the Manpack receiver. 

Overall, the data using a five-satellite constellation 
were as good or better than the data using only four 
satellites. The results for 12 May using five satellites 
(fig. 7.59) show a mean distance of 11.4 meters and a nedian 
of 10.6 meters. This compares favorably with the results 
from sets 1and2 (fig. 7.39 and 7.40), which had nean 
distances of 10.1 and 11.7 meters, respectively. The five 
satellite data from 13 May showed a very large variance due 
to the large NAVSTAR 1 URE. Therefore, data were not 
compared to remaining 13 May data. Data on 15 May for the 
five satellite constellation had a bimodal distribution 
because NAVSTAR 1 was uploaded. Once the uploading 
occurred, the peak distance was 137 meters (fig. 7.60) which 
compares favorably with 121, 241, 154, and 172 meters for 
sets 8, 9 and 10, 11, 12 and 13. Thus with a five-satellite 
constellation, positioning results were as good or better 


than the four-satellite constellation. 
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Figure 7.58. 
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Pigures 7.61 to 7.74 show the difference by fix between 
the latitude and longitude plots of the theodolite ship 
position, the differential GPS position, and the uncorrected 
GPS position. During most of the testing, the longitude of 
both methods of GPS positioning was in closer agreement with 
the theodolite ship position than was the latitude. Shown 
on figure 7.9 and 7.10 are the satellite elevation and 
azimuth data for the testing period. This data indicates 
that the satellite constellation had a larger effect on 
latitude than longitude because the predominant ranging was 
computed from one quadrant. 

The large differences between the GPS position relative 
to the theodolite Ship position were easily identified 
during the shipboard testing period. Latitude and longitude 
Blots (fig. 7.65, 7.66, 7.69 through 7.748) for sets 4, 5, 9, 
10, and 12 show the deviations of the GPS receiver positions 
from the predicted positions that were obvious to the GPS 
Operators during the shipboard testing. For sets 4 and 5, 
the stationary GPS receiver positions began to deviate 
greatly from the geoceiver position, resulting in a large 
differential correction to the GPS ship position. These 


large differential corrections caused the differential GPS 
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ViETI. CONCLOD STOWS 

The objective of the study was to determine whether the 
application of GPS in a translocated or differential mode 
provides sufficient position accuracy for near-shore 
hydrographic surveys. Several assumptions were made prior 
to the test: (1) two GPS receivers required for the 
differential mode performed the computations similarly, so 
that in the worst case a constant bias would be displayed if 
both receivers were colocated; (2) the process of updating 
the ephemeris information would not contribute a significant 
positional error to GPS in the differential mode; (3) 
reasonable accuracy could be obtained by using the output 
fron the receiver control display unLt, that isi, 
sophisticated processing methods were not required; and (4) 
the Manpack receiver results would not differ when used in 
stationary or low-dynamic environments (shipboard at 4 
knots). 

None of these assumptions proved to be valid. It was 
originally expected that the two receivers would produce 
Similar positions if both sets were operating at the same 
time over the same location. One limitation of these 


receivers, however, is that current ephemeris information is 


ae | 








not received unless the receiver has been activated prior to 
the update. {Ref. 66] For a period of time after the 
receiver is activated, the portion of the satellite message 
that contains the updated ephemeris is not acknowledged. 
Therefore, position degrades rapidly when a NAVSTAR is 
updated and the current ephemeris information is not 
utilized by the receiver. For this study, unfortunately, 
the initial times of activation of the two receivers were 
not recorded. Hence, the effect on position caused by this 
receiver parameter was not defined. 

The URE plots illustrate that the ephemeris updates did 
contribute Significantly to the range error and, 
conseguently, the resulting position error. The position 
determined by the Manpack receivers was not affected as 
seriously by the clock error in WNAVSTAR 1 when a 
five-satellite constellation was used. The inaccuracy of 
pseudo-ranges received from this satellite required frequent 
updates from Vandenburg, especially on 13 May. A receiver 
that had the option of totally eliminating the erroneous 
pseudo-range would have given the best results for ‘the 
study. Since the receivers operated independently, and the 
ephemeris information was not concurrent, the differential 


mode was effective for only a few cases. 


Lia? 








A simple correction was applied to the values output on 
the control display unit of the Manpack receiver. This did 
not permit application of tropospheric or refraction error 
corrections, other than that applied by the receiver during 
computation of the position. A more sophisticated 
prcecessing method, resulting in greater accuracy, would 
require acquisition of the pseudo-ranges and subsequent 
postprocessing to determine the receiver positions. 

The best results during the onboard portion of the study 
were obtained when the ship was allowed to drift at 
approximately 1 knot on 12 May. 

In conclusion, use of two Phase I Manpack receivers in 
the differential mode meets the reguired accuracy of 10 
meters only if the ship speed is less than 1 knot and other 
parameters such as ephemeris information and the number of 
satellites used for the position computation introduce a 
minimum amount of error. Further testing using differential 
GPS should incorporate more sophisticated processing nethods 
and more advanced receivers to achieve the survey accuracy 


reguired for nearshore hydrographic surveys. 
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APPENDIX Aw DESCRIPRION OF ERROR THEORY 





From application of error theory in the evaluation of 
positional information, it is possible t¢o establish a 
meaningful accuracy Statement subject £0 unifora 
interpretation. To provide a logical and acceptable basis 
POR COMPUTATION AND COMPARISON, POSITIONAL ERRORS ARE 
assumed to follow a normal distribution. The assumption is 
valid because positional error components generally follow a 
normal distribution pattern when sufficient data is 
available. 

The statistical treatment of errors utilizes measurable 
errors obtained in the sources of positioning information. 
Analysis of the linear components provides a two-dimensional 
expression of the accuracy of the positioning system in the 
form of an error ellipse. The use of an error ellipse is 
Palen cated by the problem of axes orientation and 
propagation of elliptical errors. Therefore, the ellipse is 
commonly replaced by circular form which is easier to use 
and understand. The linear standard errors are combined and 


converted to a circular distribution and the final accuracy 


statement is expressed in terms of circular errors. 
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The circular form assumes the circular probability 
distribution function. This function expresses the 
probability that the radial error will be egual to or less 
than the radius of the probability circle. The dispersion 
of errors within the distribution is measured in terms of 
precision indexes, each of which represents the error which 
is unlikely to be exceeded for a given probability. The 
preferred circular precision indexes are the circular 
standard error (0), the circular probable error (CPE or 
CEP), the circular map accuracy standard (CMAS), and the 
Circular near-certainty error, or 3.5 Sigma (3.50 ). For 
each precision index, the given probability is 39.35%, 50%, 
90%, and 99.78%, respectively. 

Each of the precision indexes can be converted to one of 
the others by using the following table: 


Circular Error Conversion Factors 


39. 35% 99. 8% 





39.35% | 1.0000 lest ay & 2.1460 3.5000 


50% 0.8493 1.0000 158227 2.9726 


90% 0.4660 0.5486 1.0000 1.6309 


99.78% | 0.2857 0.3364 0.6131 1.0000 
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The rapid approximation for determining the circular 

standard error is 

CSE =. 0.5000 (09x + Oy) 

where ° x and Sy are linear standard errors obtained from the 
sources of positioning information, and therefore are the 
axes of the error ellipse. For a tuuly (cazccular 
distribution, where %c is identical to the circular standard 
error, o x andoy are equal andthe angle of intersection 
between the two is 90. In all other cases, a normal 
circular error distribution is substituted for the 
elliptical distribution. The substitution is satisfactory 
tor an analysis within specifiedomin/omax ratios, where°o min 
is the minor or smaller linear standard error of the two. 
Distortion in the error distribution between elliptical and 
Circular forms at lower Omin/ oO max ratios limits the use of 
the circular concept to ratios greater than 0.2. 

The CMAS is the precision index used as the J.S. 
National Map Accuracy Standard. It 1s interpreted as 
liniting the size of error which 90% of the positions will 
not exceed. CMAS is computed from circular standard error 
as follows: 


CMAS = 2.14600Cc 
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The two-dimensional SEkor daistributions for positions 
determined by the ARTEMIS system and the theodolite network 
are described by CMAS. The accuracy of the specific 
positions is expressed by a statement of probability and 
error magnitude. The accuracy statement does not mean that 
the error in position is exactly the value shown in figure 
5.1, rather it expresses the probability that the error in 


position will not be larger than the error given. 
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